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SCHAEFFER’S REGULARITY THEOREM FOR SCALAR CONSERVATION 
LAWS DOES NOT EXTEND TO SYSTEMS 

LAURA CARAVENNA AND LAURA V. SPINOLO 


Abstract. Schaeffer’s regularity theorem for scalar conservation laws can be loosely speaking formu¬ 
lated as follows. Assume that the flux is uniformly convex, then for a generic smooth initial datum the 
admissible solution is smooth outside a locally finite number of curves in the (i, x ) plane. Here the term 
“generic” is to be interpreted in a suitable sense, related to the Baire Category Theorem. Whereas other 
regularity results valid for scalar conservation laws with convex fluxes have been extended to systems 
of conservation laws with genuinely nonlinear characteristic fields, in this work we exhibit an explicit 
counterexample which rules out the possibility of extending Schaeffer’s Theorem. The analysis relies on 
careful interaction estimates and uses fine properties of the wave front-tracking approximation. 
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1. Introduction 


We are concerned with mild regularity properties for systems of conservation laws in one space 
dimension, namely equations in the form 


d t u + d x [G(u)\ =o. 


(i.i) 

In the previous expression, the unknown U attains values in ' 

The flux function G : M. N -A 
1973, Schaeffer |24j established a regularity result (see Theorem 0 below) that applies to scalar 
conservation laws. This paper aims at showing that this result does not extend to the case of 
systems, i.e. to the case when IV > 1. 


and depends on (t, x ) G [0, +oo[xR. 
is of class C 2 . If N = 1, we call (O) scalar conservation law In 


When N > 1, system (1.1) is called strictly hyperbolic if the Jacobian matrix DG(U ) admits N 
real and distinct eigenvalues 

(1.2) X^U) < ••• < X N (U). 

We term r\(U), ..., r/v(t7) the corresponding right eigenvectors of DG(U ) and we say that the 7-th 
characteristic field is genuinely nonlinear if 

(1.3) VXi(U) -niU) > c> 0, for every U G M. N 

and for some suitable constant c > 0. In the previous expression, • denotes the standard scalar 
product in 1^. If the left hand side of (1.3) is identically zero, then the 7-th characteristic field is 
termed linearly degenerate. 


In the present paper we deal with the Cauchy problem posed by coupling (1.1) with the initial 
datum 


(1.4) U(0, ■) = Uq 

and we refer to the books by Dafermos m and Serre 
systems of conservation laws. In particular, it is well-known that, even if Uq is smooth and (1.1) is 


for a comprehensive introduction to 


a scalar conservation law, the classical solution of (1.1), (1.4) breaks down in finite time owing to 


the formation of discontinuities. The Cauchy problem (1.1), (1.4) can be interpreted in the sense of 


distributions, but in general distributional solutions fail to be unique. In the attempt at restoring 
uniqueness, various admissibility conditions have been introduced: we refer again to HUES ] for an 
overview. 

In the following we briefly go over some well-posedness and regularity results for systems of 
conservation laws. We firstly focus on the scalar case N = 1. The celebrated work by Kruzkov |l9j 
establishes global existence and uniqueness results in the class of so-called entropy admissible solu¬ 
tions of the Cauchy problem (1.1), (1.4) under the assumption that Uq G L°°. Regularity properties 


of entropy admissible solutions have been investigated in several papers: here we only mention some 
of the main contributions and we refer to mm for a more complete discussion. First, a famous re¬ 
sult by Oleinik [[22] establishes the following regularization effect: when the flux G G C 2 is uniformly 
convex, for every t > 0 the solution U(t, •) has bounded total variation, namely U(t, •) G BV(M), 
even if Uq is only in L°°. More recently, Ambrosio and De Lellis [2] improved Oleinik’s result show¬ 
ing that, except at most countably many times, the solution U(t, •) is actually a special function of 
bounded variation, namely U{t, •) G SBVfR)] we refer to [3[ § 4] for the definition of SBV(M). This 
is a regularizing effect of the nonlinearity. A result due to Schaeffer [[24] , moreover, states that for 
a generic smooth initial datum the admissible solution of the Cauchy problem is even better than 
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this: it develops at most a locally finite number of discontinuity curves, see Theorem o below. 
In the following statement, we denote by <S(R) the Schwartz space of rapidly decreasing functions, 
endowed with the standard topology (see [23] p. 133] for the precise definition). 


Theorem 1.1 (Schaeffer). Assume that N = 1 and that the flux G is smooth and uniformly convex, 
namely G"(JJ ) > c > 0 for some constant c > 0 and for every U E R. 

Then there is a set $ C <S(R) that enjoys the following properties: 
i) ^ is of the first category in <S(R), namely 


(1.5) 


5 = |^J Ck, Ck is closed and has empty interior, for every k. 


k =l 


ii) For every Uq E 5(R)\3', the entropy admissible solution of the Cauchy problem (1.1), (1.4) 
enjoys the following regularity. For every open bounded set ft C [0, +oo[xR there is a finite 
number of Lipschitz continuous curves Ti,..., T m C R 2 such that 


u ec°°(n \U 


The curves I'd,... ,T m are usually termed shocks. We briefly comment the above result. First, 
the assumption that G is uniformly convex can be relaxed, see for instance Dafermos |12|. See 
also m for recent related results. Second, a characterization of the set 5 can be found in a paper 
by Tadmor and Tassa [ 26] , Third, the result is sharp in the sense that one cannot hope that the 
regularity holds for every smooth initial datum. More precisely, even in the case G(U ) = U 2 /2 
several authors constructed initial data in <S(R) that develop infinitely many shocks on compact 
sets; see for instance the counter-example exhibited by Schaeffer himself [24, § 5]. Among recent 
works, we mention the construction by Adimurthi, Ghoshal and Veerappa Gowda |T]. 

The present paper aims at discussing whether or not Schaeffer’s Theorem 1 1.1 1 extends to systems, 
i.e. to the case when N > 1. Investigating whether or not the number of shocks is (generically) finite 
is motivated not only by intrinsic interest, but also by applications. In particular, knowing that the 
limit solution admits at most finitely many shocks simplifies the study of several approximation 
schemes. As an example, we recall that the proof of the convergence of the vanishing viscosity 
approximation in the case when the limit solution has finitely many, non interacting shocks was 
provided by Goodman and Xin m and it is considerably simpler than the proof in the general 
case, which is due to Bianchini and Bressan [6]. 

We now recall some well-posedness and regularity results for systems of conservation laws. The 
pioneering work by Glimm m established existence of a global in time, distributional solutions 
of the Cauchy problem (1.1), (1.4) under the assumptions that the system is strictly hyperbolic, 
that each characteristic held is either genuinely nonlinear or linearly degenerate and that the 
total variation of the initial datum Uq is sufficiently small. Uniqueness results were obtained in 
a series of papers by Bressan and several collaborators: we refer to the book [8] for an overview. 
In the following, we call the solution constructed by Glimm admissible solution of the Cauchy 
problem (1.1), (1.4). Note that this solution can be also recovered as the limit of a wave front¬ 
tracking approximation |8j and of a second order approximation [6], 

Several regularity results that apply to scalar conservation laws with convex fluxes have been 
extended to systems of conservation laws where every vector held is genuinely nonlinear (i.e. con¬ 
dition (1.3) holds for every i = 1 See, for instance, the works by Glimm and Lax [16] . 

Liu [21] and Bressan and Colombo |9J for possible extensions of the decay estimate by Oleinik 
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Moreover, the SBV regularity result by Ambrosio and De Lellis |2] has been extended to the case 
of systems, see Dafermos |13] for self-similar solutions, Ancona and Nguyen [4] for Temple systems 
and Bianchini and Caravenna [7j for general systems where every characteristic held is genuinely 
nonlinear. 

The main result of the present paper states that, contrary to the results by Oleinik |22j and 
Ambrosio and De Lellis [2], Schaeffer’s Theorem 0 does not extend to systems. 


Theorem 1.2. There are a flux function G : 
53 C 5(M) that enjoy the following properties: 


a compact set K C [0,+oo[xM and a set 


system (1.1) is strictly hyperbolic and every characteristic field is genuinely nonlinear, 


namely (1.2) holds and moreover property (1.3) is satisfied for every i = 1,2,3. 


ii) The set 03 is non empty and open in 5(M). 

iii) For every Uq £ 03 the admissible solution of the Cauchy problem (1.1), (1.4) has infinitely 


many shocks in the compact set I\ 

Some remarks are in order: 

• in the statement of the above theorem by shock we mean a Lipschitz continuous curve 
x = T(t) at which U is discontinuous. 


The Baire Theorem implies that any set of the first category (1.5) has empty interior. Since 


the set of “bad data” 03 is open and non empty, it cannot be of the first category and hence 
Theorem [L2] provides a counter-example to the possibility of extending Schaeffer’s Theorem 
to the case of systems. 

By looking at the explicit construction one can infer that 03 satisfies the following further 
requirement. For every Uq € 03, the total variation of Uq is sufficiently small to apply the 
existence and uniqueness results in mm • This means that the counter-example provided 
by Theorem |1.2| belongs to the same class where we have well-posedness. 

Our construction is explicit, in the sense that we provide an explicit formula for the function 
G, the compact K and the set 03, see (2.3), (5.51) and the construction in § [4] 


• Our construction shows, as a byproduct, that a wave-pattern containing infinitely many 
shocks can be robust with respect to suitable perturbations of the initial data. 

• Our counter-example requires 3 dimensions, namely N = 3. It is known that 2x2 systems 
are usually much better behaved than higher dimension systems, see for instance the dis¬ 
cussion in O § XI], An interesting question that is still to be addressed is whether or not 
Scheffer’s Theorem extends to (suitable classes of) 2x2 systems [^] 

To conclude, we briefly outline the proof of Theorem |1.2| The set 53 will be basically obtained by 
considering small IT 1 ’ 00 perturbations of a certain function U. The main point in the proof is then 
constructing G and U in such a way that 


i) when U$ = U the admissible solution of the Cauchy problem (1.1)-(1.4) develops infinitely 
many shocks, and 

ii) the same happens when Uq is a small perturbation of U. 

We choose as flux function G a particular representative of a family of fluxes introduced by Baiti and 
Jenssen [5]. Note that in |5J the authors exhibit a wave-pattern containing infinitely many shocks. 
Actually, the original wave-pattern in J5J contains large amplitude waves, but the construction can 


We thank Alberto Bressan for this remark. 
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be adapted to obtain a wave-pattern with small total variation. Although we use several results 
established in [5], our analysis is quite different from the one in (5). More precisely, there are three 
main challenges in adapting the construction in [5j for our goals : 

• we need to show that the wave-pattern in [5] can be exhibited by a solution with smooth 
initial datum: this issue is tackled by relying on the notion of compression wave, see § 4.2.1 


A much more severe obstruction is the fact that the wave-pattern in [5] is a priori not robust 
with respect to perturbations. We refer to the discussion at the beginning of § 4.3 for a 


more detailed explanation. Here we just point out that, owing to this lack of robustness, 
we have to introduce a more complicated construction than the original one in |5]. Even 
in the case when the initial datum is exactly U, the structure of the admissible solution is 
much more complex than the one considered in f5]. 

The analysis in [5] relies on the construction of explicit solutions. In our case, computing 
explicit solutions is prohibitive and hence we argue by introducing a wave front-tracking 
approximation. We perform careful interaction estimates to gain precise information on the 
structure of the approximate solution and we eventually pass to the limit by using fine prop¬ 
erties of the wave front-tracking approximation established by Bressan and LeFloch |10] . 


The paper is organized as follows. In § [2] for the reader’s convenience we go over some previous 


results. More precisely, in § 2.1 we recall some of the main properties of the wave front-tracking 


approximation, while in § 2.2 we introduce the Baiti-Jenssen system and recall some of the main 
properties. In 


system. 


3] we establish preliminary estimates on admissible solutions of the Baiti-Jenssen 
we construct the function U. 


In 


we establish the proof of Theorem 


1.2 


In particular, we show that the solution of the Cauchy problem with initial datum U develops 
infinitely many shocks and that this behavior is robust with respect to perturbations of U. 

For the reader’s convenience, we collect the notation of this paper at Page [44] 


2. Overview of previous results 

For the reader’s convenience, in this section we go over some previous results that we will need 
in the following. More precisely, we proceed as follows: 

we quickly summarize the wave front-tracking algorithm |8] and we fix some notation, 
we introduce the Baiti-Jenssen system and we discuss some of its properties. 



2.1. The wave front-tracking approximation algorithm. In this paragraph we briefly go 
over the version of the wave front-tracking algorithm discussed in [8] (see in particular Chapter 7 in 
there). We refer to [8] and to the books by Dafermos 112 § 14.13] and by Holden and Risebro [18j 
for a more extended discussion and for a comprehensive list of references. Also, in the following 
discussion we assume that each characteristic held is genuinely nonlinear (i.e., that (1.3) holds true) 
because this hypothesis is satisfied by our system. 

We first introduce some notation. We recall that the i-wave fan curve through U is 


( 2 . 1 ) 


Di[s,U) 


.Rj[s, U] s > 0 
5i[s,C7] s < 0. 
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In the previous expression, Ri is the integral curve of r) passing through U, namely the solution of 
the Cauchy problem 


( 2 . 2 ) 


dRj -> (TT\ 

*T = ri(U) ' 


Ri[0,U] = U. 


Also, we denote by Si the z-Hugoniot locus, i.e. the set of states that can be joined to U by a 
shock of the z-family, namely by a z-shock. The speed of the shock can be computed by using the 
Rankine-Hugoniot conditions. We call the absolute value |s| strength of the shock between U and 
$[«,£/]• 

We are now ready to outline the the construction of the wave front-tracking approximation. We 
fix a small parameter v > 0 and we denote by U v the wave front-tracking approximation. The 
final goal is to show that when v —> 0 + the family U v converges to the admissible solution of the 
Cauchy problem (1.1)-(1.4). The main steps to construct U u are the following (we refer to [8j § 7] 
for a detailed discussion): 


i) we construct Uq , a piecewise constant approximation of the initial datum Uq. 

ii) At each discontinuity of Uq we solve the Riemann problem between the left and the right 
state by relying on the Lax Theorem f20]• We want to define U v in such a way that U u (t, •) 
is piecewise constant for almost every t > 0. Hence, we replace the rarefaction waves 
in the Lax solution of the Riemann problem with a suitably defined piecewise constant 
approximation. The resulting approximate solution is called accurate Riemann solver. 

iii) We repeat the above procedure at each discontinuity point of Uq and we define U v by 
juxtaposing the approximate solution of each Riemann problem. In this way, U u is piecewise 
constant and has a finite number of discontinuity lines. By a slight abuse of notation, we 
call rarefaction waves the discontinuity lines corresponding to rarefactions. We can also 
introduce a notion of strength for the rarefaction wave (see 0 Chapter 7] for the technical 
details). 

iv) Let us consider the point at which two waves (i.e., discontinuity lines) interact (i.e. cross 
each other). The interaction determines a new Riemann problem, which is solved by using 
the same procedure as in step 2. above. In this way we can extend the wave front-tracking 
approximation U v after the first interaction occurs. 

v) In principle, it may happen that the number of discontinuity lines of U v blows up in finite 
time: this would prevent us from defining U u globally in time. The number of discontinuities 
can blow up if for instance U v contains a wave pattern like the one illustrated in Figure [3] 

vi) To prevent the number of discontinuities from blowing up, we introduce the so-called non 
physical waves. The exact definition is quite technical and it is given in |]8) §7.2], but the 
basic idea is the following. We introduce a threshold p, v and we consider an interaction 
point. If the product between the strengths of the incoming waves is bigger than /j, v , then 
we use the accurate Riemann solver defined at step 2. If it is smaller, we use a so-called 
simplified Riemann solver. The simplified Riemann solver involves a minimum number of 
outgoing waves. Basically, all the new waves are packed together in a single non physical 
wave, which travels at a faster speed than any other wave. 

vii) The analysis in [8j §7] shows that, by relying on a suitable choice of the approximate and of 
the simplified Riemann solver, of the approximate initial datum Uq and of the threshold //„, 
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one can prove that the approximate wave front-tracking solutions U u converge as u —> 0 + 
to the unique admissible solution of the Cauchy problem (1.1), (1.4). 


2.2. The Baiti-Jenssen system. In this paragraph we recall some results from j5]. More pre¬ 
cisely, we proceed as follows. 


|2.2.1| we introduce the explicit expression of the Baiti-Jenssen system and we comment on it. 

we recall the explicit expression of the eigenvalues and we go over the structure of the wave 
fan curves. 


2 . 2.2 


2.2.1. The system. We introduce the Baiti-Jenssen system. We fix p G]0,1[ and we define the 


function F n 


(2.3) 


by setting 


F V (U) := 


^ 4[(v — l)u — w] + rjpi(U) 

v 2 

^4 ju(u - 2 )u - (v - 1)tu j + ppziU) ^ 


In the above expression, u, v and w denote the components of U, namely U = ( u,v,w ). The 
functions p\ and p% are given by 

1 


(2.4) Pi{U) = - j[tc — ( v ~ 2) u ] 2 — [w — mi] 2 ] | = 2 uw — 2 u 2 (v — 1), 

(2.5) Pz{U) = ^ ju[u> — (v — 2 )u] 2 — (v — 2)[w — uu] 2 j = w 2 — u 2 (v — 2)v. 
In the following we are concerned with the system of conservation laws 

(2.6) d t U + d x [F v (U)\= 0, 


which we term Baiti-Jenssen system. Two remarks are here in order. First, (2.3) is exactly system 
(3.11) in |5j provided that we choose e = rj, g{y) = v 2 , a(v) = v, b{y) = v — 2, c = 4. The reason 
why we only consider a particular representative of the class of systems considered in [5] is because 


we want to simplify the analysis and the exposition. Indeed, some parts of the proof of Theorem 1.2 


are already fairly technical and we have decided to keep the rest as simple as possible. However, 
we are confident that our argument can be extended to much more general classes of systems. 

Second, the celebrated existence and uniqueness results mm mentioned in the introduction 
imply that there are constants C > 0 and <5 > 0 such that, if Uq is a compactly supported function 
satisfying 

TotVar Uq < S, 


then the Cauchy problem obtained by coupling (2.6) with the condition 17(0, •) = Uq has a unique, 
global in time admissible solution which satisfies 

TotVar U(t, ■) < C TotVar Uq , for every t > 0. 


In principle, both 5 and C depend on p. However, by looking at the proof of the convergence of 
the wave front-tracking approximation one realizes that C and S only depend on bounds on F ^ and 
its derivatives of various orders. Since all these functions are uniformly bounded in r/, then we can 
choose C and 6 in such a way that they do not depend on p. In the following, we will let i] vary 
but we will always assume that the function U attains values in the unit ball, namely \U\ < 1. 
This will be a posteriori justified because we will choose a compactly supported initial datum with 
sufficiently small total variation. 
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2.2.2. Eigenvalues and wave fan curves. We now recall some features of system (2.3) and we refer 
to 0 pp. 841-843] for the proof. First, the eigenvalues of the Jacobian matrix DF V (U ) are 

(2.7) Ai(C7) = 2rj[w — (v — 2)u\ — 4 < A 2 (17) = 2u < A 3 (17) = 2r) [w — vu\ + 4. 

Note that 

(2.8) — 6 < Ai(?7) < —^ < — 2 < A 2 (U) < 2 < 3 < \ 3 (U) <5 if \U\ < 1 and 0 < 77 < — 
and hence in particular 


(2.9) 


|Ai(£/)|, |A 2 (?7)|, |A 3 (C/)| <6 if \U\ < 1 and 0 < V < -. 


Note that (2.8) implies that the system is strictly hyperbolic if \U\ < 1 and 0 < 77 < 1/4. Note 
furthermore that 2 is a Lipschitz constant of each eigenvalue if \U\ < 1 and 0 < r] < 1/4. The first 
and the third right eigenvectors are 


( 2 . 10 ) 


n{U) = 


and r* 3 (t/) = 



respectively. The explicit expression of the second eigenvector is not relevant here. Note however 
that the assumption of genuine nonlinearity is satisfied since 


( 2 . 11 a) 

and 


VAi (U) ■ fi (U) = 4? ? > 0, VA 2 (J7) • 7 ^ 2 (17) = 2 > 0 


( 2 . 11 b) 


VA 3 (C7) • r 3 (U) = -4 ?? < 0 . 


Note that (2.11b) implies (1.3) provided that we change the orientation of r 3 . Owing to (2.10), 
the 1- and the 3-wave fan curve through U = (u, v , w) are straight lines in the planes v = v. More 
precisely, 


( 2 . 12 a) 


( 2 . 12 b) 


D^a-U] 


Ds[r-U} 


( 

y(y 


( (y + 
v 

\va + 




= U + an{U) 


= U + rr 3 (U) 


U + ar± (v ), 


U + rr 3 ( v). 


Owing to (2.11), we have that 


• if a < 0, then U and Di[cr;(u,v,w)] are connected by a 1 -shock. If a > 0, then U and 
-Di[er; (u,v,w)\ are connected by a 1-rarefaction wave. 

• if t < 0, then U and D 3 [a] (u,v,w)] are connected by a 3 -rarefaction wave. If r > 0, then 
U and D 3 [a\ (u,v,w)] are connected by a 3- shock. 


To understand the structure of the second wave fan curve through U we use the following simple 
observation, which for future reference we state as a lemma. 
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Lemma 2.1. Assume that U = ( u,v,w ) is an admissible solution of the system of conservation 


laws (2.6). Then the second component v is an entropy admissible solution of the scalar conservation 
law 


(2.13) 


d t v + d x [v 2 ] = 0. 


Proof. Lemma 2.1 was used in [5], but we provide the proof for the sake of completeness. Owing 


to the analysis in [6j (see in particular Theorem 1 and § 15 in there), the admissible solution U can 
be recovered as the unique limit e —> 0 + of the second order approximation 

d t U e + d x [F v (u e )} =ed 2 xx u £ . 

We then conclude by considering the second component and recalling that the entropy admissible 
solution of a scalar conservation law is the unique limit of the vanishing viscosity approximation 
(see [TH § 6.3]). □ 


By combining (2.11) with Lemma 2.1 and by recalling that the flux in (2.13) is convex we 


conclude tha we can choose the parametrization of D 2 in such a way that 

• if s < 0, then U and D 2 [s,U] are connected by a 2-shock and the second component of 
D 2 [s, U] is v + s < v 

• if s > 0, then U and U] are connected by a 2-rarefaction wave and the second com¬ 
ponent of Z ?2 [s 5 U\ is v + s > v. 

3. Preliminary results concerning the Baiti-Jenssen system 



when the left and the right states satisfy suitable structural assumptions. 

Both parts will be used in § [5] in the analysis of the wave-front tracking approximation of a general 
class of Cauchy problems. 

3.1. Analysis of 1-3 interactions. In this paragraph we consider the interaction between a shock 
of the first family, i.e. a 1-shock, and a 3-shock. More precisely, we term Ui, U m and U r the left, 
middle and right state before the interaction, respectively (see Figure [I] left part). In other words, 

(3.1) U m = D 3 [t, U e ], U r = D 1 [a, U m \ 


for some r > 0, cr < 0, where D 1 [•] and D 3 [-] are given in (2.12). 

We now want to solve the Riemann problem between U( (on the left) and U r (on the right). We 
recall that the 1- and the 3-wave fan curves are just straight lines in planes where the v component 


is constant, see (2.12). The slope of the lines only depends on v. This implies that the 1- and the 


3-wave fan curves commute and the solution of the Riemann problem between Ui (on the left) and 
U r (on the right) contains no 2-wave. In other words, the following holds. We denote by U' m the 
middle state after the interaction (see again Figure [I] left part). From (3.1) we get 

(3.2) U' m = D\[a, U e ], U r = D 3 [t, U'J. 
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FIGURE 1. A 1-3 interaction (left) and a 1-1 interaction (right). The value of v is 
constant across each interaction 


3.2. Analysis of 1-1 and 3-3 interactions. Owing to the particular structure (2.12) of the 1 


and 3-wave fan curves, the incoming shocks in 1-1 interactions and 3-3 interactions simply merge. 
In particular, no new wave is produced. More precisely, we have the following: we focus on 3-3 
interactions and we refer to Figure [l] for a representation. We term Ug, U m and U r the left, middle 
and right state before the interaction, respectively. In other words, 

U m — O 3 [Tt. Ug], U r — D 3 ]r r , U m ] 


for some rg, T r > 0. Owing to ( 2. 12| ) , we have U r = D\ [rg + T r , Ug] and hence the only outgoing 
wave is a 3-wave. The analysis of 1-1 interactions is completely analogous. 

3.3. Analysis of 1-2 and 2-3 interactions. In this paragraph we explicitly discuss the interaction 
of a 1-shock with a 2-shock. The analysis of the interaction of a 2-shock with a 3-shock is completely 
analogous. Lemma [3.1| below can be loosely speaking formulated as follows: if 77 and the strength of 
the incoming shocks are sufficiently small, then the outgoing waves are three shocks (and hence, in 
particular, no outgoing wave is a rarefaction). Also, we have a bound from below and from above 
on the strength of the outgoing shocks. Note that a result similar to Lemma |3.1| is established 


in [5j: the novelty of Lemma 3.1 is that we have a more precise estimate on the strength of the 


outgoing 3-shock, compare the left part of (3.6) with [5, eq. (5.9)]. Also, in the case of Lemma 3.1 


we restrict to data with small total variation. The proof of Lemma 3.1 is provided in m and is 


based on perturbation argument: one firstly establishes Lemma 3.1 in the case when 77 = 0 and 
then considers the case 77 > 0 . 

To give the formal statement of Lemma 3.1 we introduce some notation. We term Ug, U m and 


U r the left, middle and right state before the interaction, respectively. See Figure [2] left part, for 
a representation. In other words, 

(3.3) U m = D'i[s, Ug], U r = D\[(t, U m \ for some s < 0 , er < 0 . 

Also, we denote by U^ and U!^ the new intermediate states after the interaction, namely 

(3.4) U' m = D\[a', Ug], 17" = D 2 [s', U' m ], U r = D 3 [t, t/"] 


for some a', s' and rGl. Here is the formal statement of our result. 
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Figure 2. A 1-2 interaction (left) and a 2-2 interaction (right). 


Lemma 3.1. Assume that (3.3) and (3.4) hold. Then 

(3.5) s' = s. 

Also, there is £ > 0 such that the following holds. If \Ug\, |s|, |oj < 1/4 and 0 < r/ < e, then 

(3.6) — 2|cr| < a' < — ^ and ——su < r < lOscr. 

3.4. Analysis of 2-2 interactions. In this paragraph we recall a result from m concerning the 
interaction between two 2-shocks. As usual, we term Uf>, U m and U r the left, middle and right state 
before the interaction. We refer to Figure [2] right part, for a representation. Lemma 3.2 can be 
loosely speaking formulated as follows. Fix a constant a > 0 and assume that Ui, U m and U r are 
all sufficiently close to some state (a, 0, —a). Then the outgoing waves are three shocks. 


The proof of Lemma 3.2 is given in im and it is divided into two parts: we firstly establish the 
result in the case rj = 0 by relying on the explicit expression of the 2-wave fan curve. We then 
extend it to the case r] > 0 by using a perturbation argument. 

Here is the formal statement. 

Lemma 3.2. There is a sufficiently small constant e > 0 such that the following holds. Fix a 
constant a such that 0 < a < 1/2 and set U^ := (a, 0, —a). Assume that 

\Ue — t/“| < sa, s±,S 2 < 0, |si|, |S 2 1 < ea, 0 < g < ea . 

Assume furthermore that 


U r = D 2 

Then there are a < 0 and r > 0 such that 


s 2 , D 2 [si, Ui] 


(3.7) 


U r = D 3 


r, D 2 [si + s 2 ,Di[a, Ui]] 


3.5. The Riemann problem with well-prepared data. In this paragraph we discuss the struc¬ 


ture of the solution of Riemann problems with “well-prepared” data. More precisely, Lemmas 3.3 


and 3.4 below state that, under suitable structural assumptions on the constant states U~ and U + , 
the solution of the Riemann problem is obtained by juxtaposing three shocks (and hence, in par¬ 
ticular, it contains no rarefaction wave). In 5 5.3 we will use these results to discuss the wave-front 


tracking approximation of the initial datum for a general class of Cauchy problems. 
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Lemma 3.3. There is 0 < e < 1 such that the following holds. Fix Uj E M 3 such that \Uj\ < 1/2. 
Let T] i, r 2 j and r 3 j be the vectors 

(3.8) ru := ri(f7j), r 2 j := r 2 (Ui), r 2I := r 3 (C//). 

If U~,U + E M 3 safe/y 

(3.9) |tr-E/j|<e 

and 


(3.10) |C/ + — Z7 _ + hri j + 6r 2 / — 6r 2 /| < eb 
for some 0 < b < e, then the following holds. There are r, a and s such that 

(3.11) 0 < r < 26, -26 < a < 0, -26 < s < 0 

and 

(3.12) U+ = D 3 [r,D 2 [s,D 1 [a,U~}\ 

Proof. First, we point out that, if e is sufficiently small, then ( |3.10 ) implies that 

(3.13) U + - U~ = -biru - 6 2 r 2 / + 6 3 r 3 / 
for some 6i, 6 2 , 6 3 satisfying 

(3.14) 


1 3 

-6 < &i, 6 2 , b 3 < -6. 


Next, we use the Local Invertibility Theorem and we determine r, s and a satisfying (3.12). Owing 
to the regularity of the inverse map, we can infer from (3.13) and (3.14) that 

(3.15) |cr| + |s| + |t| < Cb. 

Here and in the rest of the proof, C denotes a universal constant. The precise value of C can vary 
from line to line. Next, we recall that the wave fan curve D 1 satisfies (2.12) and we introduce the 
notation 

(3.16) U' m = £>i[<7, U~ } = U~ + ari(U~) = U~ + ar u + alf^TT) - f u 
Also, we term 

(3.17) t/" := D 2 [s, U' m \ = U' m + sr 2I + s \r 2 (U' m ) - r 21 1 + \d 2 [s, U' m } — U' m — sr 2 {U' m ) 


By using (3.12) and the explicit expression of the wave fan curve D 3 (see (2.12)) we arrive at 
(3.18) 

U + = U~ + aru + sr 2I + rr 3I 

hs (to - ^3 1 


+ a 


n(U ) — fu + s r 2 (U' m ) - r 2I + D 2 [s,U' m ] - U' m - sr 2 (U' m ) 


+ T 


IZ(a,s,T,U ) 


We recall that r 2 (U^) is the derivative dD 2 [s,U^ n \/ds computed at s = 0. By using (3.9), (3.15), 
we obtain that the rest term TZ can be controlled as follows: 

| TZ(a, s, r, U~)\ < Cb(e + e) + Cb 2 + Cbe 


(3.19) 


< Cbe. 
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To establish the last inequality, we have use the assumption that 6 < e. Next, we compare (3.18) 
with (3.13) and by using (3.19) we deduce that 

\h + cr\ + \b 2 + s| + |£>3 — r\ < Ceb. 


Owing to (3.14), this implies (3.11) provided that e is sufficiently small. The proof of the lemma 
is complete. □ 


We only sketch the proof of the following lemma because it is similar to Lemma 3.3 Note, 
furthermore, that Lemma 3.3 can be recovered from Lemma 3.4 by taking the limit £ —> 0 + . 
However, we decided the give the complete statement and proof of Lemma |3.3| to highlight the 


basic ideas underpinning Lemmas |3.4| and 3.5 


Lemma 3.4. There is 0 < e < 1 such that the following holds. Fix Uj E M 3 such that \Uj\ < 1/2. 
Let r 1 1 , r 2 / and r 3 j be the same vectors as in (3.8 ). Assume that Uj , U~ , V~, U + E M 3 , and b, £ E M 
satisfy the following conditions: formula (3.9) holds and moreover 

(3.20) \V~ — U~\ < 0<6<e, 0<£< Veb. 

Assume furthermore that 


£’ 


(3.21a) 

either 

\u + - u 

(3.21b) 

or 

\u + -u 

(3.21c) 

or 

\u + — u 


Then (3.12) holds for some r, a, s such that 

(3.22a) 0 < r < 26, —26 < s < 0, —26 — £ < a < — £ 

(3.22b) 0 < r < 26, —26 — £ < s < —£, -26 < a < 0 

(3.22c) £ < t < £ + 26, -26 < s < 0, -26 < a < 0 


if (3.21a) holds 


if (3.21b) holds 


if (3.21c) holds. 


Proof. We only consider the case when (3.22b) holds since the analysis of the other cases is analo¬ 
gous, but simpler. We first rewrite (3.21b) as 

(3.23) \U + -U~+br 1 (U-) + (b + ^r 2 (U-)-br 3 (U-) + n 1 (f,U I ,U-,V~) <6/4, 

where the term TZ\ is defined by setting 


(3.24) 


£, £7/, U~, V~) :=b r 2I -r 2 {U~) +6 r u - ri(£7“) 


D 2 [-^V-]-V~ A ir 2 (V~) 


+ 6 r 3I - r 3 (U ) 

+ dr 2 (V~)-r 2 (U-) 


Owing to (3.9) and (|3.20|), it satisfies 
(3.25) 


\lZi(b, £, £7/, U~, H“)| < Ceb + C£ 2 + Cf\V~ - U~\ < C^b 

Here and in the rest of the proof, C denotes a universal constant. Its precise value can vary from 
line to line. Next, we use the Local Invertibility Theorem to determine t, s and a satisfying (3.12). 
Owing to the regularity of the inverse map, we have 


|cr| + |s| + |r| < C(b + £). 


(3.26) 
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We define U' and U!^ as (3.16) and (3.17) and by arguing as in the proof of Lemma 3.3 we conclude 


that (3.12) implies 

U + = U~ + afi(U~) + sr 2 (U~) + rr^(U~) 


(3.27) 


+ s 


+ D 2 [s,U' m }-U' m -sr 2 {U' m ) + r r 3 (C/") - f 3 (U~) 


7Z 2 (a,s,r,U ) 


By using (3.26) we obtain 
(3.28) 


\7Z 2 (a,s,T,U )| < C{b + £f 
Finally, we compare (3.23) and ( |3.27 ) and we use (3.25) and (3.28) and we obtain 
(3.29) \a + b\ + \s + b + £\ + |r - b\ < b/4 + Cj^b + C{b + £f • 


By using the inequality £ 2 < eb, we eventually arrive at (3.22b) 


□ 


By arguing as in the proof of Lemma [3.4[ we establish the following result. Note that Lemmas 3.3 


and|3.4|can be both recovered as particular cases of Lemma 3.5 


Lemma 3.5. There is 0 < e < 1 such that the following holds. Let m be a Borel probability measure 


on M. Fix Uj E R 3 such that \Uj\ < 7/2. Let ru, r 2 j and r^j be the same vectors as in (3.8). Fix 
U ~, U + E R 3 and assume that 

(3.30) \U~ — Uf < e. 

Assume, furthermore, that the the functions 


V~ 


b,i,f 2) ^3 : R —>• [0, +oo[. 


satisfy the following conditions for m-a.e. z E R: 

(3.31a) 

(3.31b) 

Finally, set 


0 < b(z) < e, 0 < ii(z) < y sb(z) for i = 1,2, 3, 
li(z) + 6(2) + 6W 1 I V~(z) - U~ | < Jeb(z). 


(3.31c) b= b(z)dm(z), fi= £i(z)dm(z), 6 = / £ 2 (z)dm(z), 6 = / £, 3 (z)dm(z). 

JR Jr JR JR 

and assume that 

U+-U~- I l D :i 


(3.31d) 


j {D 3 [Uz),D 2 [-i 2 (z),D 1 [-i 1 (z),V-(z)}}] -V~(z)}dm(z) 


+ br\ j + br 2I - br 3I 


< b/4. 


Then (3.12) holds for some r, a, s such that 
(3.32) — 26 — £i < ct < —£i, 


-2 b-£ 2 <s< -6, 


£.3 <T < i 3 + 2b. 
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4. Construction of the counter-example 


In this section we start the construction of the set of “bad data” 53 as in the statement of 
In other words, we want to construct 23 in such a way that i) 03 is open in the <S(M) 


Theorem 1.2 


topology and ii) for every initial datum in 53 the solution of the Cauchy problem develops infinitely 
many shocks in a compact set. Loosely speaking, we will construct 53 as a ball (in a functional 
space) centered at a particular function. What we actually do in this section is hence to construct 
an initial datum U satisfying the following requirements: first, the solution of the Cauchy problem 
with initial datum U develops infinitely many shocks. Second, this behavior is robust with respect 
to sufficiently small perturbations in the Sobolev space VF loo (M). As we will see in 
key step to establish Theorem 


1.2 


5.8 


this is the 


4.1 


4.2 


To construct U we proceed according to the following steps, 
we go over the construction of a wave pattern with infinitely many shocks. This construction 
is basically the same as in [5]. 

we show that this wave pattern can be obtained from a Lipschitz continuous initial datum. 
However, this does not conclude the construction of U. Indeed, at the beginning of § 4.3 


we explain that in principle it it may happen that, if we take a very small perturbation of 
the initial datum, the solution of the Cauchy problem does no more develop infinitely many 
shocks. In other words, the wave pattern constructed in 
respect to perturbations. 


4.1 and S 14.21 is not robust with 


4.3 


we modify the construction given in § 4.1 and in § 4.2 in order to make it robust with respect 
to perturbations. We eventually obtain an initial datum U and Proposition 


4.4 


states that 

the solution of the Cauchy problem with initial datum U develops infinitely many shocks 
and that this behavior is robust with respect to sufficiently small W loc perturbations. The 
proof of Proposition 4.4 is provided in § [5] 


In the rest of the present section we always assume that the parameter g in (2.3) is sufficiently 
small to have that Lemma |3. 1 1 applies. 


4.1. A wave pattern with infinitely many shocks. In this paragraph we exhibit a wave pattern 
containing infinitely many shocks. The construction is basically the same as in [5], however we recall 
it for the reader’s convenience . 

Lemma 4.1. Fix q > 0 and assume that Uj, Un, Uni £ K 3 satisfy the following properties: 

i) the solution of the Riemann problem between Uj (on the left) and Un (on the right) contains 
3 shocks and the strength of each shock is smaller than 1/4. 

ii) The solution of the Riemann problem between Un (on the left) and Um (° n the right) 
contains 3 shocks and the strength of each shock is smaller than 1/4. 

iii) The following chain of inequalities holds true: vj > vjj > vui- 


Then the admissible solution of the Cauchy problem obtained by coupling (2.6) with the initial datum 


(4.1) 


W(x) := < 


Uj x < — q 
Un —q<x<q 
U in x > q. 


contains infinitely many shocks. 

We refer to Figure [3] for a representation of the wave pattern contained in the solution of the 


Cauchy problem obtained by coupling (2.6) with the initial datum W 


















16 


L. CARAVENNA AND L. V. SPINOLO 



-q 


q 


Figure 3. The solution of the Cauchy problem obtained by coupling system (2.6) 
with the initial datum (4.1) 


Proof. First, we observe that, owing to property (iii) in the statement of Lemma 4.1 
(4.2) speed 2 [Uj , U n \ = vj + v n > v n + v in = speed 2 [i 7 / 7 , U m \ 

In the previous expression, we denote by speed 2 [C//, Ujj\ the speed of the 2-shock in the solution of 
the Riemann problem between Uj (on the left) and C /77 (on the right). In other words, the 2-shock 
that is generated at the point (t, x ) = (0, —q) is faster than the 2 -shock that is created at the point 
(t,x) = (0 ,q) (see Figure[3]). 

Next, we observe that the first interaction that occurs is the interaction between the 3-shock 
generated at x = — q and the 1-shock generated at x = q. see again Figure [3] Owing to the 
analysis in § 3.1 those two shocks essentially cross each other and, most importantly, no 2-wave 


is generated. After this interaction, the 1-shock generated at x = q interacts with the 2-shock 
generated at x = —q. Owing to Lemma 3.1 this interaction produces three outgoing shocks and 
the speed of the outgoing 2-shock is the same as the speed of the incoming 2-shock, which is the left 
hand side of (4.2). Also, the new 1-shock generated at this interaction will hit at some later time 


the left 2-shock: owing to Lemma |3.1[ this interaction produces three outgoing shocks. The new 
3-shock will then interact with the right 2-shock, producing three outgoing shocks. This mechanism 
is repeated infinitely many times between t = 0 and the time t = T at which the 2-shocks generated 
at x = —q and x = q interact, namely 

(4.3) T =———. 

vj - v m 

Note that, in general, owing to the nonlinearity, it may also happen that for instance two 3-shocks 
interact at some point on the right of the right 2-shock. However, owing to 
simply merge and no 2- or 3-waves are generated. 


3.2 these two shocks 
□ 


4.2. Shock creation analysis. This paragraph aims at showing that the wave pattern in Figure[3] 
can be exhibited by a solution starting from a Lipschitz continuous initial datum. More precisely, 
we establish the following result. 
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Lemma 4.2. There is a sufficiently small constant e > 0 such that the following holds. Fix q = 20, 
and Uj G M 3 such that \Ui\ < 1/2. Let wfR satisfy 0 < ui < £ and let Un and Uni be the states 
defined as follows: 


(4.4a) 

and 

(4.4b) 


U n ■= Ih 


Uni ■= D 3 


ui,D 2 [- ui,Di[-uj,Ui]] 


u,D 2 [-u},Di[-u,Un]\ 


Then the states Uj, Un and Uni satisfy the hypotheses of Lemma |/.1[ Also, there is a Lipschitz 


continuous initial datum such that the solution U of the Cauchy problem obtained by coupling (2.6) 
with this initial datum satisfies 7/(1, x) = W(x), where W is the same as in (4.1). 


The fact that the states Uj, Un and Uni satisfy the hypotheses of Lemma 4.1 follows from the 
remarks after formula (2.12), so we are left to prove the second part of the lemma. The proof is 
organized as follows. Since we will use the notion of compression waves in § 4.2.1 we briefly go 
over this notion for the reader’s convenience. In § 4.2.2 we give a technical lemma. In § 4.2.3 we 


eventually complete the proof of Lemma 4.2 


4.2.1. Compression waves. Consider a general, strictly hyperbolic system of conservation laws (1.1). 
We term Ri[s, U] the integral curve of r t passing through U, i.e. the solution of the Cauchy prob¬ 
lem (2.2). Assume that the i-th characteristic held is genuinely nonlinear, say VA i(U) ■ TjfU) > 0 
for every U. Let U_ := U\ for some negative s < 0 and observe that the function 

U x< \{U)t 

Ri[s, Uj x = X t (Ri[s, U])t , s < s < 0 
U x > Xi(U)t 

is a smooth solution of the conservation law on ] — oo, 0[xM and at t = 0 it attains the values 


(4.5) 


Ucw {ft %) - < 


U(0,x) = 


U x < 0 
U x > 0 


We term the function U cw defined as in (4.5) a compression wave. Loosely speaking, compression 
waves can be regarded as the backward in time analogous of rarefaction waves. 

4.2.2. A technical lemma. First, we make a remark concerning the structure of the integral curves 
R\, R 2 and R 3 of system (2.3). Owing to ( 2.10| ), we have the equalities 

(4.6) R 3 [a, U] = D x [a, U ), R 3 [r, U] = D 3 [r, U]. 


The proof of Lemma 4.2 is based on the following result. 

Lemma 4.3. Assume that the hypotheses of Lemma J^.2 are satisfied and that Un and Ujji are 


defined by (4.4). If the constants in the statement of Lemma\f.2\is sufficiently small, then 

cr, R 2 [s,D 3 [t,Ui]} 


(4.7a) 


U n ■.= D i 


for some r > -uj, s < — uj and a < — uj. Also, 

2 2 2 


(4.7b) 


Uni ■= D 1 


a*,R 2 [s*,D 3 [T*,U n }] 


for some t* > -uj, s* < - uj and a* < — uj. 

2 2 2 





























18 


L. CARAVENNA AND L. V. SPINOLO 


Proof. We only give the proof of (4.7a), since the proof of (4.7b) is entirely analogous. 

We basically proceed as in the proof of Lemma 3.3 First, we point out that (4.4) implies that 
|Ui — Uu | < Cuj < Ce. Here and in the rest of the proof C denotes some universal constant. Its 
precise value can vary from line to line. 

By using the Local Invertibility Theorem, we infer that the values of r, s and a are uniquely 
determined by imposing (4.7a). Also, we have 

(4.8) |t| + |oj + |s| < Cuj. 

We are left to prove that r > -uj, s < — -w an< ^ a < — —w. We introduce some notation: we define 
the states U' and Uf by setting 

U'-.= Di[u,Ui\, u! := D 3 [t,Uj]. 


By using (4.4) we infer 

U n = D 2 [-uj, U'] + l or 3 (v n ) 

= U'~ ujr 2 (U') + {d 2 [-uj, U'} -U' + ojr 2 ([/')} + ujr 3 (v u ) 

( 4 . 9 ) = Ui - uri(vi) - ojr 2 (U') + | D 2 [-u, U’] - U' + uf 2 (U')} + ur 3 (v n ) 

= Uj — ojri(vj) - uor 2 (Ui ) + uor 3 (vi) 

+ u{r 2 (Ui) - r 2 (t/')} + | D 2 [-u, U'] - U' + ur 2 (U')} + io^r 3 (v n ) - r 3 (u/)|. 


Note that 
(4.10) 


to 


r 2 (U /) - r 2 {U') + D 2 [—uj, U'] - U' + ur 2 (U') + u f 3 (v n ) - r 3 (u/) 


< CJ 1 . 


By using ( 4.7a| ) and by arguing as before we obtain 
Uu = Ui + rr 3 {vi) + sr 2 (Ui) + crfi(vi) 

+ s{r 2 (Uf) - i r 2 (H/)| + |i? 2 [s, Uf] - Uf - sr 2 (t/')} + 0 {f\ (pii) - n 
where, owing to (4.8), 

(4.12) s[r 2 {Uf) - r 2 (C/j)} + R 2 [s, Uf] - Uf — sr 2 {Uf) + cr|ri(u//) — ri(u/)| <Cu 2 . 

By comparing (4.9) and (4.11) and recalling (4.10) and ( 4.12| ) we obtain that 

|t — cu|T|sT£u|t|u'Tcc , | a Coj 2 . 

Since uj > 0, this implies that r > —uj, s < an d G < — -uj provided that e (and hence uj) is 
sufficiently small. □ 


4.2.3. Proof of Lemma f.2 We are now ready to complete the proof of Lemma 4.2 

We fix uj > 0 and Uj 6 M 3 , \Ui\ < 1/2. We term Un and Uni the states satisfying (4.4). 
We determine the values a, s, r, a*, s *, r* by using (4.7a) and (4.7b), respectively. Owing to 
Lemma 4.3 we have that a < 0, s < 0 and r > 0 and hence we can define the function U(t,x ) by 
“juxtaposing” six compression waves like (4.5). More precisely, we introduce the following notation: 

(4.13) U! := D 3 [r,U r ], U" := R 2 [s,Uf), U* := D 3 [r*,U n ], U** := R 2 [s*,U*} 
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For t E [0,1) we define the function U(t, x ) by setting 


(4.14) U(t,x):=< 


U r 





X 

< 

—q + a 3 (c//) • (t- 

1) 


DUj] 

if there is 0 < q < 

t: 


X 

= 

—q + A 3 (D 3 [^, Ur}) 

• it 

-1) 

m 

-q + A 3 GZO 

■ (t~ 

1) 

< 

X 

< 

-q + A 2 (H / ) • (t - 

1) 


R^M 1 ] 

if there is s < < 

0: 


X 

= 


• it- 

-1) 

Z7" 

-q + \ 2 (U") 

■ it- 

1) 

< 

X 

< 

-q + X 1 (U! , )-(t- 

1) 


D^U"} 

if there is a < < 

0: 


X 

= 

-q + XriD^U"}] 

>•(* 

-1) 

U n 

—q + Ai (Un) 

■ it- 

1) 

< 

X 

< 

q + XaiUnW-l) 



D&,U n ] 

if there is 0 < < 

t*: 


X 

= 

q + As(T>3^, Un]) ■ 

it~ 

1) 

U* 

q + HH*) 

■ it~ 

1) 

< 

X 

< 

q + X 2 (U-) ■ (t - 1) 



R2[^U*] 

if there is s* 

< q < 

: 0: 


X 

= 

q + A 2 (7?2 [?,!/*])• 

it~ 

1) 


q + HU**) 

• it- 

1) 

< 

X 

< 

q + Ai (U**) - (t - 1) 


Dik,U**] 

if there is a* 

< ? < 0: 


X 

= 

g + Ai(Diktr]) 

■ it~ 

-1) 

Um 





X 

> 

q + Ai (Um)) ■ it - 

1) 



Note that the above function is well defined because 


Ai (U„) A 3 (Un)(t - 1) + q. 


Indeed, q = 20 > 12 by assumption and \Xi(Un)\, | A 3 (C/ 77 ) | < 6 owing to (2.9). 


Note furthermore that U(t,x ) is a locally Lipsc.hitz continuous function on [0,l[xl and that 


17(1, x) = W(x), where W is the same functions as in (4.1). This concludes the proof of the lemma. 


4.3. A more robust initial datum. We firstly introduce our analysis with some heuristics. The 
analysis in the previous paragraph shows that if the initial datum is given by the same smooth 
function 1/(0, •) as in ( |4.14 ), then the solution of the Cauchy problem exhibits a wave pattern like 
the one in Figure [3] and hence, in particular, develops infinitely many shocks. However, the above 
behavior is not robust with respect to perturbations of 1/(0, •). The main obstruction that might 
prevent the formation of infinitely many shocks is the following. We recall that the strength of the 
shocks generated at time t = 1 at the points x = q and x = —q is small, more precisely it is of the 


order w < 1. By applying the second interaction estimate in (3.6), we conclude that the strength 
of the 1- and 3-shocks bounced back and forth between the two 2-shocks is weaker and weaker 
as one approaches the intersection point between the two 2 -shocks, i.e. the tip of the triangle in 
Figure [3] This means that, no matter how small a perturbation wave is, if it hits the triangle at a 
point sufficiently close to the tip it might happen that the perturbation is bigger than the shocks 
it meets. This might prevent the formation of infinitely many shocks because it might happen that 
the perturbation annihilates the shock it meets. 

In order to make the initial datum more robust with respect to perturbations we add to 17(0, •) 
the function T defined in § 4.3.2[ which is monotone in the direction of the eigenvectors. Very 
loosely speaking, the heuristic idea underpinning this construction is that in this way only shocks 
come into play, and no rarefactions. This is made rigorous in § [5] by considering the wave-front 
tracking approximation of the solution: we prove that the presence of the function T implies that 
at t = 0 the wave-front tracking approximation contains only shock waves. This will be the first 
step in the analysis that will allow us to conclude that the solution of the Cauchy problem develops 
infinitely many shocks, and that this behavior is robust with respect to perturbations. 
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We are left to make one last remark: by looking at the explicit expression of U we realize that 
there are three compression waves that interact at the point (t,x) = (1, —q) and other three that 
interact at the point (t,x) = (l,q). In § 4.3.1 we modify the datum U( 0, •) by distancing the 


compression waves one from the other. Loosely speaking, this will imply that the corresponding 
shocks will form at time t = 1 and then they will interact at some later time. This will simplify the 
perturbation analysis because it will rule out the possibility that the compression waves interact 
with each other before the corresponding shocks have formed. 

This paragraph is organized as follows: 

we modify 17(0, •) by distancing the compression waves one from the other. 


4.3.1 


4.3.2 we construct the function T “monotone in the direction of the eigenvectors” 


4.3.3 we eventually define the initial datum U in such a way that the solution of the Cauchy 
problem develops infinitely many shocks and that this behavior is robust with respect to 


perturbations. See Proposition 4.4 


4.3.1. Compression waves separation: definition ofV. We firstly introduce some notation. We fix 
a sufficiently large p > 0 (its precise value will be discussed in the following, see (4.20e)), we recall 
that the parameter q = 20 is the same as in the statement of Lemmas |4.1| and 4.2 and we set 


(4.15a) q:=g + 3, p := q — 3. 

We also introduce the following notation: 


^:=]-p,-q- A 3 ([//)[, 

fJV t :=]-q-X 3 (U!),-q-X2{U!)[ 

(4.15b) 94 m :=] — p — A i(U n ),p- A 3 (U n )[, 

%:=}p-\ 3 (U*),q-\ 2 (U*){ 

:= ]q — a 2 (IZ**),q — Ai(C/**)[, 
94 r :=]q “ Ai {Um),p[ 


94f := ] ~ q — A 3 (?7/), —q — A 3 (C/ 7 )[, 

:=] — p — MQZ"), -p - M(Un)[, 
94? := ]P — A 3 (I7ff),p — A 3 (?/*)[, 

94 ? :=]q-X 2 {U*),q-X 2 (U**)[, 

K Hq-AiOT*),q-Ai(E/zff)[, 


We also define the open sets 91 c and 9\ w by setting 
(4.15c) 94 c = 9\ e U % U m" U U U U 94" 

and 

(4.15d) 9\ w = 9lf U 9lf U 94^ U 94^ U 94^ U 94,., 


respectively. To give an heuristic interpretation of the above notation we point out that, if we had 
q = p = q, then the intervals in (4.15) would be the same as in the right hand side of (4.14). In 
particular, we would have that the function U{ 0, •) is constant on 94 c and has a nonzero derivative 


on 94 w . 

To construct the function V, we fix the parameters 5 and uo and we set 


(4.16) 


Ur := (<S,0,-<S). 
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ftt 


ft} ft!} 


ft} 


fti 


ftr 


ft} 

fti. ft" 


ft! 


ftr 


ft w 

ftr 



-P-Ai (U n ) 



P 

q — Xi(Ujji) 

q-Ai (£/**) 
q-HU**) 
q — A 2 (i7*) 

P - A 3 (tn 


FIGURE 4. Intervals defined in Equations (4.15) 


We determine the values cr, s, r, cr*, s*, r* by using (4.7a) and (4.7b), respectively. Finally, we 
determine U 1 , U". U_* and U_** by using (4.13). We now define the function V :] — p, p[—> M 3 in 


such a way that V is a 3-compression wave on ft} U ft}., a 2-compression wave on ft} U ft}, and a 
1-compression wave on ft} U 93*. More precisely, we set 


(4.17) 


V(x) := < 


U r 





x £ ftp 


if 

there 

is 

0 

u 

V 

O’ 

V 

uf 





x £ ft} 

RzfaU'] 

if 

there 

is 

s 

o 

V 

Lr 

V 

u" 





x £ ft'} 

D^U"] 

if 

there 

is 

a 

A 

A 

o 

U n 





X £ ftjn 

U n ] 

if 

there 

is 

0 

< g < t* : 

U* 





X £ ft' r 

R2[q,U*j 

if 

there 

is 


' < < 0: 






x £ ft} 

D^U**] 

if 

there 

is 

cr 

* <g <0: 

Uni 





x £ 9v r 


Note that if we had q = p = q, then V would coincide with the function f/(0, •) defined as in (4.14). 


4.3.2. Monotonicity in the direction of the eigenvalues: definition of 4U We fix the parameters 
Cc > 0 and ( w > 0 and we define the function 41 :] — p, p [—> M 3 by requiring that 4r(0) =0 and that 


(4.18) 


-Qnr - ( c r 2 i + ( c r3i if x £ ft c 
-C wfu - ( W V 2 i + (wT3i if x e ft w 


In the previous expression, we used the notation r\j = fi(Uj), r- 2 i = r 2 (Uj) and r* 3 / = V 3 {Ui). 


4.3.3. Definition of the initial datum U. We now define the Lipschitz continuous function U : M —> M 3 
by setting 


(4.19) 


U(x) 


4> _ (x) x < —p 

< V(x) + 4i(x) —p<x<p 
4> + (x) x > p. 
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FIGURE 5. The solution of the Cauchy problem with initial datum the function V 


defined as in (4.17) 


In the above expression, the function V is as in (4.17), the function T is defined in § 4.3.2 and the 
functions <3? _ , 4> + : M —> M 3 are Lipschitz continuous and defined in such a way that the function U 
is continuous and compactly supported. We also require that each component of <3? _ (a;) and 4 > + (.t) 
is monotone. 


We can now state the main result of the present section. Proposition 4.4 below states that 


i) the solution of the Cauchy problem obtained by coupling (2.6) with the initial datum 
U(0,x ) = U has infinitely many shocks; 

ii) this behavior is robust with respect to sufficiently small perturbations of the initial datum. 
Proposition 4.4. Fix q = 20. Let 0 < e < 1 and fix the parameters 

(4.20a) 5 := e, ( w := e/2, 


(4.20b) 

(4.20c) 


V := e , 
Cc:=e 9 , 


U) := £ 


r :=e w /2, 


Note that by combining the above choices with (|4.3|) and (|4.4|) we get 
(4.20d) 

We also require 


~ 20 
T = 


(4.20e) 


6 

p := 12T + 40 = 40 ( + 1 


Consider the same function U as in (4.19). If the constant e is sufficiently small, then, for every 
initial datum Uq such that 


(4.21) 


\\U 0 ~ U\\ w i<x, < r, 
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the admissible solution of the Cauchy problem obtained by coupling system (2.6) with the initial 


datum U{ 0, •) = Uq has infinitely many shocks in the bounded set ]0, 2T[x] — 2q, 2 q[. 


The proof of Proposition 4.4 is the most technical part of the paper and it is given in § [5] The 


main result of the present paper, namely Theorem 1.2 follows as a corollary from Proposition 4.4 


see 


5.8 


5. Proof of the main results 


In this section we establish the proof of Theorem 1.2 and Proposition 4.4 More precisely, in 


5.8 we show that Theorem |1.2| follows as a corollary of Proposition 4.4 The rest of the present 


section is devoted to the proof of Proposition |4.4| Since the proof is fairly technical and articulated, 
we provide a roadmap in § |5.1[ The proof is established in the remaining paragraphs. 


5.1. Proof roadmap. In this paragraph we provide the proof outline and we discuss the basic 
ideas underpinning the analysis in the following paragraphs. 

We start with some heuristic considerations. We recall that the function V : M —> M is defined 


as in (4.17). The qualitative structure of solution of the Cauchy problem with initial datum V 


is illustrated in Figure |5} by the time t = 1, six shocks have formed. More precisely, moving 
from the left to the right there are a 3-shock, a 2-shock, a 1-shock, a large interval where the 
solution is constant and then again a 3-shock, a 2-shock and a 1-shock. These shocks interact at 
some later time and produce a wave pattern with infinitely many shocks. The initial datum Uq 
is obtained from V by adding the function T and the perturbation JJq — U, which is IT 100 small, 


see (4.19), (4.21). Loosely speaking, the goal of the following paragraphs is to show that adding T 


and Uq — U to the initial datum does not affect too much the qualitative structure of the solution of 
the Cauchy problem and, in particular, does not jeopardize the formation of infinitely many shocks. 
Since computing explicit solutions is prohibitive, we rely on the wave-front tracking approximation. 
The proof of Proposition |4.4| is organized as follows: 

We make some preliminary remarks that will be used in the following paragraphs. 

We introduce the wave front tracking approximation U v (u is the approximation parameter) 


5.2 


5.3 


of the solution of the Cauchy problem obtained by coupling the Baiti-Jenssen system (2.6) 


with the initial datum U{ 0, •) = Uq satisfying (4.21). In particular, we construct a piecewise 


constant approximation of the initial datum and we discuss the waves that are generated 
at t = 0. A feature that will be very useful in the analysis at the following paragraphs is 
that at t = 0 only shock waves are generated. This is the reason why we introduced the 


function T, monotone in the direction of the eigenvectors, see (4.18), (4.19) and the analysis 
54L21 and § [57343] 


m 


5.4 We carry on a qualitative analysis of the waves of the wave front-tracking approximation 
U v . In particular, we split the wave generated at t = 0 in two groups: group A comprises 
the waves that will contribute to the formation of six “big shocks” like in the solution of the 
Cauchy problem with initial datum V. Group B comprises all the other wave generated 
at t = 0, which in the following will be regarded as perturbation waves. In § |5.4| we also 
introduce groups of waves generated at interactions occurring at times t > 0. They will 
also be regarded as perturbation waves in the following. Note that perturbation waves are 
important, even if they are small, because they contribute to the formation of infinitely 
many shocks. 
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5.5 


We establish quantitative bounds on the total strength of the waves belonging to the various 
groups introduced in 


5.4 


5.6 


We eventually establish the results concerning the shock formation, see Lemmas 5.3 and |5.4[ 
In particular, we show that the wave front-tracking approximation U v contains six “big 
shocks” like the solution of the Cauchy problem with initial datum V, see the discussion 
in § [57672] 

In particular, we firstly provide a 


5.7 We eventually conclude the proof of Proposition 4.4 


bound from below on the number of shock fronts in the wave front-tracking approximation 
U u , see Lemma 5.5 Next, we pass to the limit v —> 0 + and we conclude that the number 
of shocks of the limit solution is infinite on a given compact set. The limit analysis relies 
on fine properties of the wave front-tracking approximation established by Bressan and 
LeFloch HOl. 


We conclude this paragraph with two technical remarks. First, as pointed out in § |2.1| in this paper 
we use the version of the wave front-tracking approximation discussed in the book by Bressan |8j. 
This version involves the use of two kinds of procedures to solve wave interactions: the accurate 
Riemann solver and the simplified Riemann solver. Whether one or the other is used depends on the 
product of the strength of the incoming waves, see the discussion at the beginning of § 5.7.1 and the 
analysis in |8j Chapter 7] for more detailed information. To simplify the exposition, in § 5.4, § 5.5 
and § 5.6 we pretend we always use the accurate Riemann solver. The fact that there are actually 


two kinds of solvers is taken into account in § 5.7 


Second, to simplify the notation in the following we denote by 0(1 ) any quantity which is 
uniformly bounded and bounded away from 0, namely there are universal constants c, C > 0 that 
satisfy 

0 <c<0( 1) < C. 


5.2. Preliminary considerations. In this paragraph we collect various remarks that we will use 


in the following. First, we fix £ > 0 sufficiently small so that Lemmas 3.1 3.2 3.3, 3.4 4.2 4.3 

5 = e, 


apply. Next, we recall formula (4.20): 


p = 0(1)uj 
_2 


-i 


T = <D(l)u 


-l 


OJ = £ 


VI - VIII = U 

C W = e/2, V = e 2 , Cc = e 9 , r = e 10 /2. 

This in particular implies 

(5.1a) ( w u + ( c p + rp <£lo < £ 3/4 ( w r] 

(5.1b) r < £^ c < £u < £( w 

We will use the above inequalities in the following. 

We recall that the intervals ..., 94 r are as in (4.15) and that the function V is defined in (4.17). 
By construction, we have 


(5.2) 
and 

(5.3) 


TotVar V < 0( l)w 

1^ - UiWco < 0(l)w, ||P || C 0 < 0(1)(<5 + u). 
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We can infer from (2.11a), (4.7a), (4.7b),(4.13),(4.15) and (4.15d) that the length of 9\ w is 0(l)co 
because the length of 94f, 94|, 94]?, 94.) is 0(l)ujr] while the length of 94jj, 94): is O(l)oj. Since 'h(O) = 0, 
from ( 4.18| ), (4.20a) and (5.1) we get that 

(5.4) H^llco < TotVar T < 0(1)( w uj + 0(l)( c p < 0(1 )euj. 

Also, we recall that each component of 4> _ and 4> + is monotone and that 4?” and 4> + both attain 
the value 0. This implies that 

||$-|| c o + ||4> + || c o < TotVar 4>" + TotVar T+ < | V(-p) + ®(-p)| + \V(p) + T(p)| 

(5.5) <|^/| + |C%|+0(l)(CcP + C«;w) 

— 0(1) [5 + U! + CcP + Cio w ] • 


By recalling (4.20), (5.2) and (5.4) we conclude that 

(5.6) TotVar U < 0(1)6. 


Owing to (4.21), we have 

(5.7) 


\\U 0 - U\\ c o + TotVar(V 0 -U)< 0(l)rp 


If Uq satishes (4.21), which means that Uq is a perturbation of U, then 

(5.8) Uq{x) = V(x) + 4t(x) + [Uq(x) — U(x)] for every x e] — p, p[. 

By the explicit expression of V and by (|2.11[) we infer that | V'(x)\ < 0( l)f? _1 for every x E] — p, p[. 


By using (4.18), (4.21) and ( |5.8[ ) we arrive at 

(5.9) \Uq(x)\ < 0(l)ri~ l for every x G] — p, p[. 


By taking into account (5.3), (5.4), (4.21) and (5.1), we infer from (4.20) and (5.8) that 
(5.10) \U 0 (x) - Uj\ < | U 0 (x) - U(x )| + (x)| + \V(x) - U t \ 

< r + 0(1)eu + 0(1)lo < 0( l)e 3 for every x E] — p, p[. 


We point out that by estimates (4.20b), (5.1a), (5.2) (5.4) and (5.7) one has the bound 
(5.11) TotVar Uq < 0(1)uj < C9(l)e 3 on]— p, p[ 


Since Ui = (5,0, —5), owing to (4.20a) we arrive at 

(5.12) |V 0 (x)| < 0( l)e, for every x G] — p, p[. 

Remark 5.1. We point out that the values attained on ] — 2 q, 2g[x]0,2T[ by the admissible solution 
of the Cauchy problem are only determined by the behavior of the initial datum on ] — p , p[. This 
follows by combining our choice (4.20e) of p with the finite propagation speed, more precisely 
with (2.9). Indeed, we have 

p — 2q > 12 T > max |Aj(t/)| • 2 T. 

|U|<1,*=1,2,3 

In the following, we will only be concerned with the behavior of the initial datum on the interval 
] — p, p[. This is justified by the previous considerations and by the fact that we are only interested 
in the behavior of the solution on ] — 2 q, 2c/[x]0, 2T[. 

5.3. Wave front-tracking approximation: initial datum. In this paragraph we discuss the 
wave-front tracking approximation of the initial datum. We recall that the intervals 94^,..., 94 r are 
defined in (4.15). 
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5.3.1. Mesh definition. We fix an approximation parameter v > 0 and a mesh size h u > 0. We 

require that h v —> 0 + when v -A 0 + . We choose Xq < x\ < ■ ■ ■ < x ^ in ] — p, p[ so that 

(5.13a) (1 — e)h v < x" +1 — x" < h v for every i = 0,..., m u — 1. 

If h u is sufficiently small, one can as well assume that the extrema of the intervals 9b, • • •, 9t r are 
all contained in the set {xq, ..., x^}. Define the wave-front tracking approximation of the initial 
datum by setting 

(5.13b) Uq(x) := Uq(x^) for x E]x", x" +1 [ and i = 0,..., rn v — 1. 

We now describe the waves generated at the grid points Xq,..., x by separately considering the 

regions 9b, ■ • ■ • 9t r . 


5.3.2. Waves generated in 9t c = 9b U 91^ U 91" U 9t m U 91) U 91" U 91 r . We only focus on the analysis 
of the interval 9b because the analysis of the other intervals is entirely similar. 

We fix x" E 9b and we consider the Riemann problem between the states 

U~ := lim Uq(x) = Uq(x\_ 1 ) (on the left), U + := lim Uq(x) = Uq{x”) (on the right), 
ad«)' a;4.(a:D + 


Claim. If (5.1a) and (5.1b) hold, then the states U ,U + ,Uj := Uj satisfy the hypotheses of 


Lemma 


3.3 with the choice b = Qfxf — x"_i) 


Proof. Hypothesis (3.9) in the statement of Lemma 3.3 follows by (5.10). Next, we focus on 


hypothesis (3.10). We use (5.8) and we recall that V is constant on each connected component of 


9t c , while VH' = ( c (—ru — r* 2 / + ? 3 l)- This implies that, if b = C,c{xf — x"_ 1 ), then 

\U+ _u_ + Wu + b r 2I _ b r 2I \ = \(U 0 - U)(xf) - {U 0 - U)(xU)\ ^ - x^r), 


which owing to (5.1b) gives inequality (3.10). 


□ 


CONCLUSION: By using Lemma |3.3[ we conclude that the only waves created in the open set 
9t c are 1-, 2- and 3-shocks. In particular, no rarefaction waves are generated. Moreover, owing 
to (3.11) the total variation of all these waves is bounded by 0(l)CcP- 


5.3.3. Waves generated in 9f u , = 9t| U 9f(? U 9t| U 91^ U 91j) U 91.). We only focus on the analysis of 
the interval 9lj) since the analysis of the other intervals is entirely similar. We fix x" E 9l| and we 
consider the Riemann problem between the states 

U~ := lim Uq(x) = t/o(x"_i) (on the left), U + := lim U{f(x) = Uq{x\) (on the right). 

Claim. Assume that U~,U + are as at the previous line and that V~ := R(x"_ 1 ). Let £ > 0 be the 
strength of the 3-shock between V~ (on the left) and V(xf) (on the right), namely 


V(xf) = D 3 [C,V~}. 


If b = Cw( x i ~ xf- 1 ) } then all the hypotheses of Lemma 3.) are satisfied. 


Proof. Hypothesis (3.9) in the statement of Lemma 3.3 follows by (5.10). Next, we point out that 


the condition 0 < b < £ is satisfied provided that v is sufficiently small. Indeed, b < 0(l)£ w h u and 
h v —> 0 + when v —> 0 + . 
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To check the other hypotheses, we first recall that G 94). By combining the explicit expression 
of V (4.17) with (2.11a) and (2.11b) we infer that the derivative of V satisfies V'{x) = 0(l)?y _1 . 
This implies that £ = 0(1) (x” — x\_-^)r]~^ and hence that 

= 0(1)Ve(wV- 

Next, we plug (5.4) and (5.7) into (5.8) and we get that the first condition in ( 3.20[ ) is satisfied: 
\V--U~\ = \^>(xU) + Uo(xi-i) - U(xt i)l < 0(l)(( w u + C cP + r) 


(5.14) 


< 0(l)s( w r] < y/e-. 


To establish the last inequality we used (5.1a). The condition £ 2 < eb is satisfied because 
£ 2 = 0(1)« - <_,) V 2 = 0( l)h 2 r?- 2 < b = 0(l)s( w h u 


provided that h u is sufficiently small. Finally, we check that (3.21c) holds. We use again (5.8) and 
we recall that 'F / = ( w (—rn — r* 2 / + r 3 j) on Dl w . By using (4.21) and (5.1b), this implies 

| U+-IT- D 3 K, V~] + V~ + br u + br 2I - br 3I \ 

= I£/(,«) - Uo(xU) - VK) + V(xU) - ^«) + 

= Wo - u)(x\) - (Uo - U)(xU)\<r« - <-i) < 

All the hypotheses of Lemma [3.4| are therefore satisfied. This concludes the proof of the claim. □ 


By applying Lemma 3.4 
Conclusion: By Lemma 


and using (3.22) we arrive at the following conclusion. 

the only waves created in the intervals Ttj?, 1H 2 9\ l f and 94) 


3.4 


are 1-, 2- and 3-shocks. In particular, no rarefaction waves are generated. 

The total variation of all the 3-shocks generated in the intervals 94^, 94)! is 0(l)u and the total 
variation of all the 1 and 2-shocks generated in the intervals 94)!, 94)! is bounded by 0(l)C, w u:r}. 

The total variation of all the 2-shocks generated in the intervals 94 2 , 94 2 is 0(l)u and the total 
variation of all the 1 and 3-shocks generated in the intervals 94 2 , 94 2 is O(l)Q w 0 J. 

The total variation of all the 1-shocks generated in the intervals 94^, 94), is 0(l)u and the total 
variation of all the 2 and 3-shocks generated in the intervals 94), 94,) is 0(l)^ w ur]. 


5.4. Wave front-tracking approximation: qualitative interaction analysis. In this para¬ 
graph we split the waves of the wave front-tracking approximation U u into several groups, that 
are defined in the following. As we will see in § 5.6 and as we pointed out in the proof roadmap 


m 


5.1 the waves of group A are the waves that will contribute to the formation of a wave pattern 


similar to the one of the solution with initial datum V (see Figure |5|. The waves of groups B can 
be heuristically speaking regarded as perturbation waves. 

We now define the groups A, B, Ci, ..., C m . In § |5.3| we discussed the waves that are generated 
at t = 0. In particular, we proved that only shocks are generated at t = 0. We split these waves 
into two groups: 


• Shocks of group A: group A comprises 

— the 3-shocks generated in the intervals 94) and 94) and their right extreme; 

— the 2-shocks generated in the intervals 94) and 94) and their right extreme; 

— the 1-shocks generated in the intervals 94) and 94) and their right extreme. 
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We fix a shock i 6 A, Let V t be its strength, which is defined as in § 2.1 Due to the 


(5.15) 


conclusions at the end of § 5.3.3 the total strength of all the shocks of group A is 0(1)cj, 
namely 

J> = 0(l)w 

A 

SHOCKS OF GROUP B: group B comprises all the shocks generated at t = 0 in the interval 
] — p, p[ which are not comprised in group A. In other words, group B comprises 

— the 1, 2 and 3-shocks generated in the open interval SK C ; 

— the 1 and 2-shocks generated in the intervals 91') and 91'j: 

— the 1 and 3-shocks generated in the intervals Dx'j and 91 j: 

— the 2 and 3-shocks generated in the intervals and 91). 

Owing to the conclusions at the end of 


5.3.2 and of § 5.3.3 the total strength of all these 


shocks can be bounded by 


(5.16) 


<O(l)(p( c + 0jQ. 


ieB 


We now want to track the evolution of the shocks of groups A and B by discussing their interactions. 
For the time being, we do not take into account the fact that in some cases we have to use a simplified 
Riemann solver (see § 2.1). We will take into account the presence of non-physical waves in § 5.7 


We separately consider the following cases: 

i) We fix two shocks, i and j, and we assume that i is either a 1 or a 3-shock and j is a 2-shock. 
Just to fix the ideas, let us assume that i is a 3-shock. In this step, we do not care whether i 
and j belong to group A or B. Let Vi and Vj be their strengths and we assume that i and j 


(5.17) 


n 


interact at some point. By Lemma 3.1 all the outgoing waves are shocks. By definition, we 
still call i the outgoing 3-shock and we still call j the outgoing 2-shock. Also, the outgoing 
i belongs to the same group (A or B) as the incoming i, and the same happens for j. We 
say that the outgoing 1-shock is the new shock which is created at the interaction. This 
new shock belongs neither to A nor to B: we define a new group Ci in the following. 

We conclude by recalling some interaction estimates: let V- and Vj be the strengths of i 
and j after the interaction. By Lemma 3.1 Vj = Vj. Also, we recall [8] formula (7.31) p. 
133], which states that 

|V'-V*| <0(1)^,, 

Also, [8j formula (7.31) p. 133] implies that the strength of the new shock generated at the 
interaction is bounded by 0(l)ViVj. 

We fix two shocks, i and j, and we assume that i is a 3-shock and j is a 1-shock. Owing 
to the analysis in § 3.1 the outgoing waves are a 3-shock and a 1-shock. By definition, 
we still call i the outgoing 3-shock j the outgoing 1-shock. We say that the outgoing i 
belongs to the same group as the incoming i. The same holds for j. Finally, we recall some 
quantitative interaction estimates: we term V* and Vj, V[ and Vj the strengths of i and j 
before and after the interaction, respectively. Owing to the analysis in § |3.1[ 


(5.18) 


V' = 


Vj = V, 
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iii) We fix two shocks i and j and we assume that they are both 2-shocks. We term V* and Vj 
their strengths and we assume that they interact at some point. We set 


a := 5, 


U* :=U T = (<5,0, -<5), 


- by . 


s 2 = -V, 


and we claim that the hypotheses of Lemma 3.2 are satisfied. The hypothesis 


|Ug — | < ea holds owing to (5.10) and to the fact that a = 5 = e (see (4.20a)). 


0 < r] < ea is satisfied because rj = e 2 owing to (4.20b). 

|si|, |s 2 1 < e are satisfied. Indeed, owing to Lemma 2.1 the maximal strength of a 2- 


shock is bounded by the total variation of the v component. The total variation of the 
v component at t = 0 satisfies TotVaruo < 0( l)e 3 by (5.11). Since the total variation 


of a scalar conservation law is a monotone non increasing function with respect to 
time [ T9] . we can conclude that |si|, |s 2 | < £■ 

Lemma |3.2| states then that the outgoing waves at the interaction point are three shocks. 

We now separately consider the following cases: 

• if i belongs to A and j belongs to B, then we term i the outgoing 2-shock and we 
prescribe that it still belongs to A. Note that the strength of i after the interaction is 
Vj = Vi + Vj. We set V' = 0, in such a way that 


(5.19) 


V' + V'=V; + V,, 


iv 


• if i and j both belong to either A or B, then we proceed as follow. Just to fix the ideas, 
assume that i is the fastest shock among the two, namely i is on the left of j before the 
interaction. By definition, we still call i the outgoing 2-shock and we prescribe that it 
belongs to the same group (A or B) as the incoming shocks. We also set Vj =0, in 
such a way that (5.19) holds. 

In both cases, we say that the outgoing 1 and 3-shock are new shocks generated at the 
interaction. Note again that these new shocks belong neither to A nor to B: they will 
belong to the group Ci defined in the following. 

We fix two shocks, i and j, and we assume that they belong to same family, which can be 
either 1 or 3. Just to fix the ideas, let us assume that they are both 3-shocks. Owing to the 
analysis in § 3.2 the only outgoing wave is a 3-shock. We separately consider the following 


cases: 


if i belongs to A and j belongs to B, then we term i the outgoing 3-shock and we 
prescribe that it still belongs to A. Note that the strength of i after the interaction is 
V' = Vi + Vj. We set V\ = 0 in such a way that (|5.19|) holds. 


if i and j both belong to either A or B, then we proceed as in case iii). Just to fix the 
ideas, assume that i is the fastest shock among the two, namely i is on the left of j 
before the interaction. By definition, we still call i the outgoing 3-shock and we say 
that it belongs to the same group as the incoming shocks. We also set Vj =0, in such 
a way that (5.19) holds. 


We explicitly stress three properties following from the analysis of cases 0 > E above. First, the 
outgoing waves are always shocks. Second, new shocks are only created when the interaction 
involves at least one 2-shock. Third, the new shocks created at the interaction are either 1- or 
3-shocks, namely no new 2-shocks are created. 
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We now focus on the new shocks created at interaction points. They can only be 1- and 3-shocks, 
since by Lemma 2.1 and by the definitions in [I]) and iii) above no new 2-shock can arise. We now 
collect them into a sequence of groups C m defined by recursion on m £ N. 

• Group Cp we term Ci the group of new shocks that are generated at the following 
interactions: 

— between a shock i and a shock j both belonging to group A. 

— between a shock i belonging to group A and a shock j belonging to group B. 

— between a shock i and a shock j both belonging to group B. 

As mentioned before, Ci only comprises 1- and 3-shocks. It follows from the analysis in 


3.1 


that if two shocks of group Ci interact then either they basically cross each other or they 
merge: to label the outgoing waves at the interaction, we proceed as in case[jv]) above. Note 
furthermore that, if a shock i £ Ci merges with a shock j £ A U B of the same family, we 
term j the outgoing shock and we set V[ = 0, Vj = Vi + Vj. Hence, the only possibility for 
the generation of new waves is the one discussed at the next item. 

Group C m +i: we term C m +i the group of new shocks that are generated at interactions 
between a shock i belonging to group C m and a 2-shock j belonging to either group A or B. 
As mentioned before, C m only comprises 1- and 3-shocks. At interactions among shocks in 
C m+ i the shocks can basically either basically cross each other or they can merge. In any 
case, no new shock is created: to label the outgoing waves at interaction points, we proceed 
as in case [iv]) above. Also, if a shock i £ C m +i merges with a shock j £ C m U A U B of the 
same family, we denote by j the outgoing shock and we set V- = 0: in this way equality 


(5.19) is satisfied. 


In this way we have classified all the shocks of the wave-front tracking approximation U u . 

5.5. Wave front-tracking approximation: quantitative interaction estimates. This para¬ 


graph aims at establishing Lemma 5.2 below. In the statement, Vi denotes as usual the strength 
of the shock i. 

Lemma 5.2. There is a constant K > 0 such that, if the constant e in the statement of Proposi¬ 
tion \4.f \ is sufficiently small, then we have the following estimates: for every t > 0 


(5.20a) 

E ^ Kuj - 

(5.20b) 

E V ^(t + ) < Kus, 

*bSB 

(5.20c) 

E v i( t+ ) < ( 2A ^) 

iECm 


m+1 


Proof. We point out that owing to (5.15) and (5.16) combined with (5.1a) we can choose K in such 
a way that we have the inequalities 

(5.21) ]T V lA (t = 0) < \ku, E V ifl (* = °) < Ik(u( w + pQ < \ku>e. 

*a£A i s e B 

The shocks of groups C m , m £ N, do not exist at t = 0, but we can adopt the notation that their 


strength is 0, in such a way that (5.20c) is formally satisfied. The proof of the lemma is based 


on the following argument: we assume that estimates (5.20a), (5.20b) and (5.20c) are satisfied for 


every t < t and we show that they are satisfied for t = t. The technical details are organized in the 
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following four steps. 

0 Step 1: we make some preliminary considerations.. We first introduce a new notation: we 
denote by D the group 

oo 

(5.22) D:=BuJC m . 

771=1 


In the above expression, the groups B and C m are as in § 5.4 Note furthermore that here an in 
the following we term groups the sets A, B, C, n , D, while we use the term family as a shorthand 
for characteristic family. 


Note that by combining all the inequalities in (5.20) and (4.20b) we get that, if e is sufficiently 
small, then 

OO 

Y Vi(t+) < Kies + Y (2-Kw) 


(5.23a) 

(5.23b) 


vm+1 < 2Kus for every t < t 


iGD 771=1 

Y v *( t+ ) - Ku + 2Kuje ^ 2Kuj 

ie AUD 


for every t < t 


Note furthermore that the quantities at the left hand side of (5.20), (5.23a) can only change at 
interaction times. 

0 Step 2: we establish the bound on Xa^sA^a- Note that the only ways X^gA^a can change 
are the following interactions: 

i) Interactions where a shock i E A with strength V,; merges with a shock j with strength Vj of 


the same family and of the same group A. In this case (5.19) ensures that yA Vj A does 
not change at this interaction. For this reason, in the following we neglect these interactions. 

ii) Interactions where a shock Ia EA with strength V-; A merges with a shock jo of the same 
family but of group D. In this case V- A = Vi A + Vj D . Each shock jo of group D may have 
at most one of these interactions: let J l ff be the subset of shocks of group D that merge 
with i\. 

iii) Interactions where a shock i a EA with strength V* A interacts with a shock j of a different 


family. In this case by the interaction estimate (5.17) one has 


V' A < v tA + o(m A Vj = V lA (l + <D(l)Vj ). 

Each shock j may interact at most once with a given shock i& of a different family. 

We recall that all shocks in group A are generated at time t = 0 and we track the evolution of a given 
shock Ia GA between time t = 0 and t = t. If the shock i a only interacts with shocks j E CD 
and with waves j±, ■ ■ ■ ,jk of different families then by |n|, [m| above one has the inequality 


v iA (t + ) < i v lA (t=o)+ y I n 1 + 




jeAuD 


Note that the last factor in the above expression does not depend on ia- Also, by using the 
inequality e x > 1 + x we get 

Yl (l + 0(l)Vj) < e° (1)E >^u»b. 

jeAuD 
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We now sum over all the shocks %a EA and be obtain 

E v u(*" + )< (e v iA (t = o)+x; E v i| IT 


*a£ A 




i A eA je fA I jeAuD 

Vi.A (* = 0) + E V i ) ex P ( C> ( 1 ) 

In 6 D / \ 

The last inequality holds because two sets and are disjoint subsets of D (keep in mind 
that we are neglecting the fact two shocks EA and i' A EA can merge). We now plug the above 




inequality into (5.21) and we recall that by assumption at time t < t estimates (5.20) hold. Owing 


to (5.23) we obtain 


(5.24) 


E 


Vi A {t + ) < ^-Kcu + 2K wej exp (0(1)Kuj) < Kui 


provided that e is sufficiently small, since owing to (|4.20bl) ui = e 3 . Note that (|5.24|) implies 


that (5.20a) holds for t = t provided that (5.20b) and (5.20c) hold for t < t and (5.21) holds at 
t = 0. 

Q Step 3: we control X^ b gb^*.b- The only ways Yhi B e&Vi B can increase are the following: 

i) If a a shock is EB merges with a shock j eAUB of the same family. In this case Vi B 


does not increase owing to (5.19). 


ii) If a shock is EB with strength V* s merges with a shock j E U me isjC m of the same family. 

iii) If a shock is EB with strength V* s interacts with a shock j of a different family. 


By repeating the same argument we used in Step 2 we conclude that, if the initial estimate (5.21) 


holds and moreover (5.20) holds for t <t, then by (5.23) 


(5.25) 


E 

*b£B 


< 


E V Jfl (f = 0) + 

B 


oo 

EE 

m= 1 ieC m 


Vi exp 


oo E 


Vi 


jeAuD 


< {\ Kuj£ + eXP - Ku£ 


provided that e is sufficiently small, due to (4.20b). Inequality ( |5.25 ) implies that (5.20b) holds for 
t = t provided that (5.20) hold for t < t and ( 5.21| ) holds at t = 0. 

Q Step 4: we conclude the proof. To control ^ igC V* we firstly recall that a shock of group 
Ci can be generated when a shock j\ belonging to either group A or B interacts with a shock 72 
belonging to either group A or B. As pointed out before, the strength of the outgoing new shock is 
bounded by 0(l)Vj 1 Vj 2 . We denote by Vf the strength of the shock j E Ci at the time when the 
shock i is generated. We then have 


E < E E °( l )VnVh + E E °wvnVn + E E °wvhV ja 

ie Cl iieAj 2 eA jiGAj 2 £B ji£Bj 2 £B 

< K 2 (u 2 + u 2 e + (we) 2 ) < 2I\ 2 oj 2 . 

Next, recall that X^gCi Vi(t) can increase not only when a new shock is generated, but also: 

i) if a shock i ECi with strength V* merges with a shock j E C m , m > 1 of the same family. 

ii) if a shock i ECi with strength Vi interacts with a shock j of a different family. 
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By arguing as in Step 2, we find that if (5.20) holds for t < t and (5.21) holds at t = 0, then 


OO 


£ 

ieCi 


Vi(t + ) < ( V i + y l Vi I exp I C ’( 1 ) 

71=2 

( 2 Ku) 3 \ 



< 


2K 2 u? + 


1 - 


2^3 ) exp(C>(l)Ahu) < (2Ku) 


provided that e is sufficiently small, since owing to (|4.20b|) u = e 3 . We have thus established (|5.20c|) 


for all t > 0 when m = 1. The case when m > 1 can be handled in an entirely similar way. This 
concludes the proof of Lemma 5.2 □ 


5.6. Wave front-tracking approximation: shock generation analysis. In this paragraph we 
finally show that the in wave-front tracking approximation one can recognize a wave pattern like 
the one of the solution of the Cauchy problem with initial datum V, see Figure [5] In particular, in 


5.6.1 we establish the generation of six “big shocks”: Lemma 5.3 establishes the formation of two 


1-shocks and two 3-shoc.ks, while Lemma |5.4| established the formation of two 2-shocks, which are 
moreover approaching. In 


5.6.2 we conclude the analysis of the wave pattern generation. 


5.6.1. Shock formation: small times. We recall that the interval 94 3 is defined by formula (4.15) 
and we establish the following lemma. 

Lemma 5.3. By the time t = 6/5, some of (or all of) the 3-shocks of group A generated at time 
t = 0 in the interval 94 3 merge into a single 3-shock with strength greater than ui^/e/2. 

Proof. We first describe the idea underpinning our argument. We term jf and j r the 3-shocks that 
are generated at t = 0 at the left and the right extrema of the interval 94 3 , respectively. In Step 3 
below we show that at t = 0 these two 3-shocks are approaching. We then track the the evolution 
of jf and j r on the time interval ]0, 6/5[ and we point out that there are only two possibilities: 

• The strength of both jf and j r remains smaller than ui^/e. In this case we show in Step 5 
below that jf and j r keep approaching and they merge by time t = 6/5. We also show that 
this implies the creation of a 3-shock with strength at least O(l)co. 

• The strength of either jf or j r surpasses ujy/e at some time t E]0, 6/5[: just to fix the ideas, 
let us assume that it is the strength of jf. In Step 6 below we show that this implies that 
the strength of the jf remains bigger than co^/e/2 on the whole interval ]t, 6/5]. 

The technical details are organized as follows. 

Q Step 1: we point out that in the time interval ] 0,6/5 [ the waves of group A generated in 94 3 


can only interact among themselves and with the waves of group D (see (5.22) for the definition of 


group D). In other words, they cannot interact with shocks of group A generated in other intervals. 


To see this, we proceed as follows. We recall definition (4.15) and that the shocks of group A 
are only generated in the intervals 94 3 , 94/, 94/, 94/, 94/ and 94/. The closest interval to 94 3 is 94/ 
and the distance between the right extreme of 94 3 and the left extreme of 94/ is 


(5.26) 


- q + q - A 2 (U 1 ) + A 3 (t/') > 3 - A 2 (U!) + A 3 (C/') > 6. 


To establish the last inequality, we used (5.12) and the explicit expression of the eigenvalues, 


see (2.7). This implies that, if the constant e in the statement of Proposition 4.4 is sufficiently 


small, then in the time interval ]0, 6/5[ the 3-shocks generated at t = 0 in 94 3 cannot interact with 
the 2-shocks generated at t = 0 in 94/. 
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Q Step 2: we focus on the time t = 0 and we introduce some notation. Let xp and x r be the right 
and left extrema of namely xp = — q — As(f7/) and x r = — q — A 3 (U_')- We recall that by the 
mesh definition discussed in § 5.3.1 xp and x r are both points of discontinuity for Uq . Next, we 
define the states Uf ( 0), [//(0), U~{ 0), J7+ ((1) by setting 

(5.27a) 

(5.27b) 


UK 0) := lim UX(x), 

XJX£ 


Up~ ( 0 ) := lim Ctf (x) 

X\Xl 


U r (0) := lim Uq(x), 

XjX r 


17+(0) := lim U%(x) 

x\.x r 


We denote by jp and j r the 3-shocks of group A generated at t = 0 at xp and x r , respectively. We 
also denote by speedy ( 0 ) and speedy ( 0 ) their speed at t = 0 . 

Q Step 3: we control from below the initial difference in speed of jp and j r . More precisely, we 
establish the following estimate: 

(5.28) speedy( 0 ) - speeded) > ^ length(lKf) = ^ (a 3 (£/j) - A 3 ( V )')) • 


To this end, we point out that owing to (5.9), 
|L/( 0 )-[/+( 0 )|<O(l)^, 


|[/-( 0 )-[/+( 0 )|<O(l)^. 

V 


The explicit expression (2.7) of A 3 implies that |VA 3 1 < 0(l)r] and hence by using the above 
inequalities we arrive at 

(5.29) jspeedy(0) - A 3 (t/ £ -(0))| < 0(l)h„, |speed jr (0) - A 3 (C/+(0))| < 0(l)h u . 

Next, we use ( |4. 21 ), (5.4), ( |5.8| ), ( |5. 13b ) and the equalities V(xp) = Ui and V(x r ) = U' to get 

(5.30) |C7(0)-E7/| <0(1) (eu + r + ^j , |t/+(0) - C/'| < 0(l)(ecn + r). 

Exploiting again the equality |VA 3 | < 0(1)?/, we get that (5.30) implies 

|A 3 (£/)T(0)) - A 3 ([//)| < 0(1) (ewry + rr] + h v ) 

|A 3 (C7+(0)) - A 3 (C/j)| < 0(1) (ecnr? + rr ,). 

By plugging the above estimate into ( |5.29 ) we arrive at 

(5.31) speedy(0) - speedy(0) > A 3 (L/) - X^U 1 ) - 0(1) (ewq + rrj + h v ) 

Next, we point out that the equality |VA 3 | < 0(l)r/ implies that 

(5.32) A 3 (0/) - A 3 (t/') = 0(l)wr/, 


because by Lemma 4.3 the parameter t in (4.13) is of order co. We eventually obtain (5.28) by 
observing that terms in the last parenthesis in (5.31) are of lower order than ojij. this follows by 
recalling (5.1b), and the fact that h u 0 when v 0. 

0 Step 4: we consider the evolution of the shocks jp and j r in the time interval ]0, 6/5[. Let £// (t) 
and U7(t), f/+ (tj and U~(t ) be the left and right state at time t of jp and j r . respectively. Note 
that the above functions are piecewise constant: to define their pointwise values, in the following 
we choose their right continuous representative. One of the following two cases must occur: 
i) we have 

(5.33) \Up-(t) ~ UjT(t)\ < ujy/s, | £//(£) — U~(t)\ < u^/e for every t €]0,6/5[. 

We handle this case in Step 5 below. 
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ii) There is t e]0, 6/5[ such that 

(5.34) either |£//(f) — Uf(t)\ > ojyfe or \U+(t) — U~(t)\ > uy/e. 

We handle this case in Step 5 below. 


Q Step 5: we conclude the proof of the lemma under the assumption that (5.33) holds. 


We recall from Step 1 that in the time interval t E]0, 6/5[ both jp and j r can either merge with 


other 3-shocks of group A or interact with 1-, 2- and 3-shocks of group D (5.22), but they cannot 


interact with other 1- or 2-shocks of group A. This implies that U^ (t) and £// (t) can only change 
owing to the interaction with some shock of group D: we recall ( 5.23a| ) and we conclude that 

(5.35) | Uf(t) - C//(0)| + | U+(t) - U+( 0)| < 0 (\) uje . 

Next, we proceed as in Step 3 and by combining (5.33) with ( |5.35 ) we conclude that 

speed ± (t) - speed jr (t) > | length(94/) = ^ (\ 3 (Ui) - A 3 (C/)) for every t E]0, 6/5[ 


(5.36) 


provided that e (and hence w, owing to (4.20b)) are sufficiently small. In the previous expression, 


we denote by speedy(i) and speedy(t) the speed of jp and j r at time t. Note that (5.36) implies 
that by the time t = 6/5 the shocks jp and j r merge. By construction, this implies that all the 
3-shocks of group A generated at t = 0 in 94/ merge by time t = 6/5. In the following, we denote 
by A/ the group of the 3-shocks of group A generated at t = 0 in 94/. We follow the same argument 
as in Step 2 of Lemma |5.2| and we use the inequality 

n (i-omjs i-o(i) x; Vj, 

jeAuD jeAuD 

which is a consequence of the elementary inequality (1 — x)(l — y) > 1 — (x + y) if x,y > 0. We 
conclude that the total strength of the shocks in A/ can bounded from below, more precisely by 
recalling (5.23a) and the analysis in § 5.3.3| we have 


i£D 


Vi) £ Vi(t = 0) > 0(l)cu. 


a \ 


Vi(t)> (i-o(i); 

ieAf 

We eventually obtain that by time t = 6/5 the shocks of group A/ merge into a single shock with 
strength 0(l)u. 


Q Step 6: we conclude the proof of the lemma under the assumption that (5.34) holds. 


First, we point out that (5.34) implies that at t = t part of the waves of group A/ have merged into 


a shock of strength uiy/e. Hence, we are left to prove that this shock “survives” with a sufficiently 
large strength up to time t = 3/2. To this end, we point out that for t > t this shock can merge 
with other 3-shocks of group A/ and hence increase its strength. Also, it can interact with other 


shocks of group D: however, by following the same argument as in Step 2 of Lemma 5.2 and by 


recalling (5.23a) the strength of the shock is bounded from below by 

0(1) ^2 Vi)uy/e > UJyfs/2, 


1 - 


i&D 


provided that e is sufficiently small. This concludes the proof of Lemma |5.3 


□ 


Note that by repeating the above proof we obtain the an analogous of Lemma 5.3 holds for 
the 3-shocks of group A generated at time t = 0 in the interval 94/ and the 1-shocks of group A 
generated at time t = 0 in the intervals 94/ and 94/.. In the case of 2-shocks we have a stronger 
result. 
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Lemma 5.4. Let T be the same constant as in (4.3). The following conclusions hold true: 

i) By the time t = 6/5, all the 2-shocks of group A generated at time t = 0 in the interval 
9 \ 2 merge into a single 2-shock J 2 having strength greater or equal than 0(l)u. The same 
holds for the 2-shocks of group A generated at time t = 0 in the interval 91 2 , let jf be the 
resulting shock. 

ii) The 2-shocks jf and jf are approaching and they merge by the time t = 2 T. 

Proof. The proof of [IJ) is organized in two steps. We only discuss the 2-shocks of group A generated 
at time t = 0 in the interval 91 2 , the argument for the 2 -shocks generated in 9l 2 is completely 
analogous. 

Q Step 1: we discuss the situation at time t = 0. We recall that, owing to Lemma 2.1 the speed 
of a 2-shock j between U~ and U + is 


(5.37) 


speed,,- = v +v + . 


Next, we fix x ”, xf x G £R 2 . We denote by j t and the 2-shocks generated at t = 0 at the points 
x = x\ and x/ +1 , respectively, and by speedy. (0) their speed at t = 0. Let v' U q , v’ v be the first 
derivative of the second components of Uq and V, respectively. By combining (4.17), (4.18), (4.21) 
and (5.8) we have 

v'u 0 ( x ) + \ = \ v u 0 ( x ) - V v\ < 0(l)Cw + 0(l)r, for every x G JK? 


By using the relations Q v = e. r < e 2 , (5.13a) and (5.37), we get that, when e is small enough, the 
last inequality brings us to 


(5.38) 


speeded) - speed ii+ 1 ( 0 ) > Q - 0( l)e^ «_i - <+i) > Q - 0{l)ej 2(1 - e)h v 




Q Step 2: we show that the shocks ji and merge in the time interval ]0, 6/5[. By the 
arbitrariness of x\ and xf +1 this establishes |ij). 

We recall that the speed of a 2-shock does not change at the interaction with a 1- or a 3-shock. 
Hence, the speed of ji and can only change when they merge with a 2-shock. Three cases can 
occur: 

• the shocks ji and ji+\ merge: this proves the claim of the present step. 

• ji merges with a 2 -shock t on the left of jp. this implies that the speed of j: increases. 

• j ^)_i merges with a 2 -shock l on the right of ji +\: this implies that the speed of ji+\ 
decreases. 


If only the last two cases occur, by recalling (5.38) we conclude that 

5 

speed j.ft) - speed ji+1 (t) > speedy (0) - speed ii+ 1 (0) > -« - < +1 ) for every t G]0,6/5[. 
This implies that by the time 6/5 the shocks ji and j,+l merge, and hence concludes the proof of 

0- , n 

0 Step 3: We are now left with establishing nil), namely proving that the shocks J 2 and J 2 
(defined as in the statement of Lemma 5.4) merge by the time 2 T. To this end, we recall the 
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9\ 3 r :=}q-X 2 (U*),q-X 2 (U**)l 


explicit expression of 91/ and 91/ (4.15): 

Kj :=)- q -X 2 (U!),- q -X2(U!% 

We also introduce the following notation: we term 

• vj (t) the second component of the left state of the 2-shock created at t = 0 at the left 
extreme of 91/, 

• v ( t ) the second component of the right state of the 2-shock created at t = 0 at the right 
extreme of 91/, 

• speed^(f) the speed of the 2-shock arising at t = 0 at the left extreme of 91/. 

The functions v~(t), v+(t) and speed r (i) are similarly dehned by considering 91/. Note that vj(t) 
and vf(t) are the left, right state and speed of j/ for t > 6/5, because by (/) all the 2-shocks 
generated at t = 0 in 91/ merge by the time t = 6/5. By using an analogous argument we prove 
that v~(t), v+(t) and speed r (f) are the left state, the right state and the speed of <//, respectively. 
By combining (4.17), ( 4.21| ), (5.4) and ( |5.8[ ) with ( |5.1a| ) and (5.13b) we infer that 


(5.39) 


{[v e ( 0 ) + n/( 0 )] - [v r ( 0 ) + v+( 0 )]} - [vi - v m \ < 0(1) (ue + K). 


Next, we point out that v ( v r vif can only vary with respect to t owing to the interactions 


with 2-shocks of group B. Owing to (5.20b), this implies that 

(5.40) | vj(t) - v7(0)| + | vf(t) - n/(0)| + \v~(t) - u“(0)| + | v+(t) - n+(0)| < O(l)co£. 


By combining ( 5.39[ ) and (5.40) we infer 

speed t (t) - speed r (t) > [vj ( t) + vj(t)\ - [v~(t) + v+(t)] 
> [vi - vni} ~ 0( l)(we + h v ). 


(5.41) 


By using (5.41) and the dehnitions (4.3), (4.15) of T and 91;, we realize that the shocks Jj and J, 


merge by time 

^ [g-HU**)]-i-g-HU!)] ^ 

sup t [speedy (t) — speed,, (t)] ~ 


2q + 0(l)u: 
vi - vni - 0(1) (euj + h v ) 
l + 0(l)u ~ 


2q + 0(l)u 
2 q/T — 0(1 )(elo + h v ) 


l-0(l)(e + TK) 

To get the last equality we have used the equalities T = ^(l )^^ 1 and q = 20. Since h u —> 0 + , 
this implies that, if u> = e 3 is sufficiently small, then jj and ,// merge by the time t = 2T. This 
concludes the proof of Lemma |5.4[ □ 

5.6.2. Shock formation: wave pattern generation. By relying on the analysis at the previous para¬ 
graph, at t = 6/5 the wave-front tracking approximation U u (t, •) contains at least six “big shocks”. 
Going from the left to the right, i.e. as x increases, we encounter: a 3-shock with strength at least 


ujy/e/2 (see Lemma 5.3), a 2-shock with strength greater or equal than 0(l)u, a 1-shock with 
strength greater or equal than ujy/e/ 2, and then again 3-shock with strength at least cjy / e/2, a 
2 -shock with strength greater or equal than O(l)u), a 1 -shock with strength greater or equal than 
w-y^/ 2 . Note that the two 2 -shocks are approaching and they meet by time t = 2T. Also, the six 
big shocks do not interact on the time interval ]0, 6/5[ because the generation regions 91/, 91/, 91/, 


91/, 91/, 91/ are sufficiently separated, see (|4.15a|) and (4.151) Besides those six “big shocks” there 


are in general other waves, which however are all shocks by the analysis in § 5.4 
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5.7. Conclusion of the proof. In this paragraph we conclude the proof of Proposition 4.4 
In § 5.7.1 we take into account the presence of non-physical waves in the wave-front tracking ap¬ 


proximation. In § 5.7.2| we establish a bound from below on the number of shocks in the wave-front 
tracking approximation. Finally, in 


5.7.3 we complete the proof of Proposition 4.4 


5.7.1. Non-physical waves. In this paragraph we take into account the presence of non-physical 
waves. We firstly recall some facts about the simplified Riemann solver and we refer to [8, §7.2] 
for a complete discussion. 

First, one chooses a threshold parameter /a u > 0. We discuss the choice of p, u later in this 
paragraph, however we point out that pL u —> 0 + as v —> 0 + . The accurate Riemann solver is used 
to solve the interaction of a wave a of strength V a with a wave (3 of strength Vp in the wave 
front-tracking approximation if the product of the strengths of the incoming waves satisfies 


(5.42) 


' Vp T Rw 


If the above condition is violated, we use the simplified Riemann solver, which is defined at (8J 
p.131] and involves the introduction of so-called non-physical waves. Non-physical waves travel at 
a speed faster than any other wave and the simplified Riemann solver is defined in such a way that 
their interaction with the other waves has a minimal effect. To simplify the exposition, here we 
do not recall all the technical details and we only discuss the properties of the simplified Riemann 
solver for the Baiti-Jenssen system (2.3) that we need in the following. These properties are either 
a direct consequences of the definition of simplified Riemann solver or can be straightforwardly 
recovered by combining the definition with the features of the Baiti-Jenssen system discussed in 

SIS3I 


i) If the incoming waves are a 1-shock and a 3-shock, then the simplified Riemann solver 
coincides with the accurate Riemann solver. 

ii) If we use the simplified Riemann solver to solve the interaction between a 2-shock and a 1- 
shock (respectively a 3-shock), then the outgoing waves are a 2-shock, a 1-shock (respectively 
a 3-shock) and a non-physical wave. The value of the v component is constant across the 
non-physical wave. 

iii) If the incoming waves are both 1-shocks then the simplified Riemann solver coincides with 
the accurate Riemann solver. The same happens if the incoming waves are both 3-shocks. 

iv) If the incoming waves are both 2-shocks, then the speed of the outgoing shock is the same 
in the simplified and in the accurate Riemann solver. Also, the value of the v component 
is constant across the outgoing non-physical wave. 

v) By combing all the above features we conclude that the strength of the v component is 
always constant across non-physical waves. 

vi) If a non-physical wave interact with a 2-shock, then the speed of the 2-shock does not 
change. 

vii) The strength of each non-physical wave is at most p, v . Also, owing to the analysis in [8] 
p. 142] we can choose fi u in such a way that the total strength of non-physical waves satisfies 


(5.43) 


total strength non physical waves < u, 


where v is our approximation parameter. Owing again to the analysis in [8] p. 142], this is 
consistent with the requirement that \x v —> 0“*" as v —> 0 + . 
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We now discuss how the presence of the non-physical waves affect the analysis at the previous 
paragraphs. First, we point out that it does not affect at all the discussion on the initial datum 
in § |5.3| because by definition the simplified Riemann solver is only used at time t > 0. Next, we 
point out that the use of the simplified Riemann solver forces the total number of waves to be finite. 
In particular, there actually are fewer waves of groups Ci,..., C m than those considered in § |5.4[ 
Lemma |5.2| does not change if we take into account the presence of non-physical waves, provided 
that we say that if a group C m is empty, then the the total strength of its waves is 0. The reason why 
Lemma |5.2| does not change is because the proof is based on interaction estimates on the strength 
of waves and by definition the interaction with a non-physical wave does not change the strength of 


a shock. Also Lemma 5.3 does not change: indeed, the proof is based on the quantitative estimates 


given by Lemma 5.2 which are still valid. The further perturbation provided by the non-physical 


waves is arbitrarily small owing to (5.43) and hence does not affect the proof. Finally, Lemma 5.4 


does not change because the proof is based on estimates that, as a matter of fact, involve only the 
second component (i.e., the component v) of the wave front- tracking approximation. Owing to 
properties §,0, [v]) and[vlj of the non-physical waves of the Baiti-Jenseen systems, non-physical 
waves have basically no effect on the v component of the wave front-tracking approximation and 


hence the proof of Lemma 5.4 is still valid if we take into account non-physical waves. 


5.7.2. A bound from below on the number of shocks. This paragraph aims at establishing Lemma 5.5 
below. In the statement, Jf and J 2 are the same as in the statement of Lemma 5.4 and we denote 


by [•] the entire part. Also, fi u is the threshold to determine whether we use the accurate or the 


simplified Riemann solver, see (5.42). 


Lemma 5.5. Fix a threshold 6 > p v /bJ 2 . In the bounded set (t, x) e] — p , p[x]0, 2T[, the wave-front 
tracking approximation U v admits at least 


(5.44) n g := 

shocks j such that strength Vj of j satisfies 

(5.45) 


lof V 2 (—) 


L> 


Proof. If there were only the six “big shocks” mentioned in § 5.6.2, then the wave pattern would 


be qualitatively like the one represented in Figure [5] To understand the impact of the other waves 
and to establish (5.44) we track the evolution of the left 3-“big shock” jf, which has strength at 


least uy/e/4: when it interacts with the left 2-“big shock” Jf. We recall that the strength of Jjt is 
0(1)uj and we use estimate (7.31) in [8j p. 133]: we conclude that after this interaction the strength 
of Jf is at least 

u -f _ o(i)v^ 2 > 

After this interaction, the shock Jf moves towards the right 2-“big shock” jf. Before interacting 
with J 2 , however, j'f can interact with 1- and 3-shocks and with 2-shocks different than jf and 
J 2 . The interaction with a 3-shock increases the strength of Jf because the shock merges with 
jf. The interaction with a 1-shock does not affect the strength of Jf. We are left to consider the 
interactions with 2-shocks different than jf and J 2 . We recall that 2-shocks are only generated at 
t = 0 and that, owing to Lemma 5.4 all the 2-shocks generated at t = 0 in and 94 2 have merged 
by the time t = 3/2 to generate J 2 and J 2 , respectively. Hence, what we are left to consider are 
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the interactions of Jf with the 2-shocks that are not generated at t = 0 in !SR? U 5H 2 . Note that 


all these 2-shocks belong to group B. We recall (5.20b) and the interaction estimate (7.31) in [HI 
p. 133] and we infer that after the interaction with all these 2-shocks the strength of j'f is at least 


(5.46) 


Owing to Lemma 


Uli 


--o ( u^Ev, 


— 


^-o( i) w V/ 2 >^. 

8 w “16 


3.1 


when jf interacts with J 2 . then a 1-shock is created: by combining (5.46) 
with the fact that the strength of J 2 is G(l)w, we get that the strength of this 1-shock is at 
least 0{l)y/£uj 2 /16. Also, this 1-shock moves towards J 2 , but before reaching jf may interact 
with 1-, 2- and 3-shocks. By arguing as before, we infer that when it reaches J 2 its strength 
is at least 0(l)i/ew 2 /32. When this 1-shock interacts with J 2 , a 3-shock with strength at least 
0(l)y/suj 5 /32 is created. We repeat this argument as long as the strength of the reflected 1- or 


3-shock j satisfies (5.45), namely we can repeat it a number ng of times, where ng satisfies 

0(1)^ (!)”*> e , 


This implies (5.44). We are left to justify the fact that we used the accurate and not the simplified 
Riemann solver. Note that, owing to the inequality 9 > p u /u) 2 and since the strengths Vj 2 of J 2 


and Vj 2 of J 2 are equal to if (5.45) holds then 


Vj -max{Vj 2 , Vj 2 | > 0(l)9uj > 0(1)— > p u 

r i (jJ 


provided that oj is sufficiently small. 


simplified Riemann solver and it concludes the proof of Lemma 5.5 


This implies that we must use the accurate and not the 

□ 


5.7.3. The limit solution has infinitely many shocks. We are eventually ready to conclude the proof 
of Proposition 4.4 First, we rely on the analysis in [HI Chapter 7] and we conclude that when 


v —> 0 + the wave front-tracking approximation U u (t,-) converges strongly in L 1 1 oc (M) to a limit 
function U(t,-) for every t > 0. Also, the function U is the admissible solution of the Cauchy 


problem obtained by coupling (2.6) with the initial datum f7(0, •) = Uq. 

We are left to prove that U admits infinitely many shocks in ] — p, p[x]0,2T[. We rely on 
Lemma 5.5 and on fine properties of the wave front-tracking approximation established in |10j 


(see also j8j §10.3] for an introductory exposition). 


More precisely, we firstly point out that the function ng defined as in (5.44) satisfies 


(5.47) 


lim ng = Too 
0 ->O+ 


since u> < 1. Next, we refer to the definition of maximal 9-shock front given in [8j p.219]: loosely 
speaking, a maximal 0-shock front is a polygonal line made by consecutive shocks of the same 
family where the strength of each shock is greater or equal than 9/2 and there is at least one shock 
having strength greater or equal than 9. 

Also, we consider the “big” 2-shocks J 2 and J 2 given by the statement of Lemma 5.4 We term 
(t*,x*) their intersection point and we remark that by construction x* G [—2q,2q], Note that by 


looking at the proof of Lemma 5.5 we realize that, if 9 > p u /ui 2 , then there are at least ng shocks 
with strength bigger or equal than 9 and that cross the part of the plane between J 2 and J 2 , 
namely they intersect the vertical line x = x* at some time t < t*. 
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We now argue inductively as follows. We fix a threshold 6 ± > 0 such that ng 1 > 1, namely there 
is at least one shock jf such that the strength of jf is at a some point greater or equal than 6 \. In 
particular, jl is a maximal #i-shock front. 

Next, we fix 62 such that ng 2 — ng 1 /2 > 1: this implies that, for very v sufficiently small, U u has 
at least a shock j'f satisfying 


(5.48) 


Q 

62 < strenght £ < 


By arguing as in j8j p. 220] we infer that when v —> 0 + the shock curves jl and converge 
uniformly (up to subsequences) to two shocks of the limit function U: we term them and 
Also, the value x* converges (up to subsequences) to some limit value x*. Note that the limit 
shocks and j'f^ both intersect the vertical line x = x* and, moreover, the strength of jl 


is 


greater or equal than 6 \ and the strength of j'f. is comprised between 62 and 0\/2. This implies 
that j^ and j'f are two distinct shock curves and, hence, the limit solution has at least 2 shocks 
with strength greater or equal than 62 . 


Owing to (5.47), we can iterate the above argument: for every natural number k , there is 6 ^ > 0 


such that the limit U has at least k distinct shocks with strength greater or equal than 6^. This 
implies that U has infinitely many shocks and concludes the proof of Proposition |4.4| 


5.8. Proof of Theorem |1.2[ This paragraph aims at establishing the proof of Theorem |1.2| Before 
entering the technical details, we make some preliminary heuristic considerations. To establish 
Theorem 1.2 we need to construct a set 03 C 5(M) that satisfies condition [nj and iii) in the 


statement of the theorem. Proposition 4.4 states that, if U is the same as in (|4. 19), then the 


admissible solution of the Cauchy problem with initial datum U develops infinitely many shocks 
and this behavior is stable with respect to W loo -perturbations. Note, however, that both U and its 
!T loo -perturbations have discontinuous first order derivatives and hence they do not belong to 5(M). 
To construct 03, we mollify U to obtain a smooth function and we consider PT loo -perturbations of 
the mollified function. 

We now provide the technical details: we first introduce the notation. We fix a convolution 
kernel 0, namely a smooth function 


(5.49) 


[0, +oo[, 


c p(x)dx = 1 , cj)(x ) = 0 if |x| > 1 . 


We fix ? > 0 and we define the mollified function : M —> M 3 by setting 
(5.50) C/ ? (x) := f U(x + sz)(j>(z)dz, 


i -1 


where U is the same function as in (4.19). Note that € 5(M) since 


is compactly supported because so it is U. 
is smooth by the classical properties of convolution. 


Theorem 1.2 is a direct corollary of the following result. 


Proposition 5.6. There is a sufficiently small constant e > 0 such that the following holds. Assume 
that q = 20 and that 5, (ci r h 0J i r am ^ P are as the statement of Proposition f.f Assume 
furthermore that g < e 2 r]( c . Let he the same function as in (5.50) and set 
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Then condition [|), [wp and Hi) in the statement of Theorem 1.2 are satisfied. 

The proof of Proposition 5.6 is divided into two parts: in § 5.8. 1| we establish a technical lemma 
which loosely speaking says that Lemma |3.5| applies to the Riemann problems obtained from the 
piecewise constant approximation of Uq E 53. In § |5.8.2 we conclude the proof of Proposition 5.6 


and hence of Theorem 11.21 

5.8.1. Analysis of the Riemann problems arising at initial time. This paragraph is devoted to the 


proof of Lemma 5.7 


We assume that the hypotheses of Proposition 5.6 are satisfied and we recall that the set 53 
is defined as in (5.51). We also recall the mesh definition in § 5.3.1 we fix v > 0, h v > 0 and 
Xq < x\ < 


< xm in ] — p, p[ in such a way that (5.13a) holds. 


Lemma 5.7. Assume that the same hypotheses as in the statement of Proposition 5.6 hold true. 
Fix x't e] — p + e, p — e[ and set 


(5.52) 

and 


U~:= lim^(.T) = t/ 0 «_ 1 ), 

xfx" 


U + := lim UX(x) = U 0 (a£) 

xixv 


V (z) := V «_! + gz) , 

where V is the same function as in (4. 17). Also, let T be the same function as in ( |4.18[ ) and 
consider the functions £i,£2)^3 : [— 1,1] —> M which are defined for every z E [—1,1] by the 
equalities 

-b(z)ru -b(z)r 2 i + b(z)r 3I = T [x\ + gz) - T (x^ + gz) 

V « + gz) = D 3 \Uz),D 2 \-Uz),D l \-^(z),V-(z) 


(5.53) 


Finally, let be the same function as in (5.49) and m the C 1 -absolutely continuous measure defined 
by setting 

(5.54) m (E) := / f>(x)dx for every C l -measurable set E. 

J E 

Then all the hypotheses of the Lemma\3.5\ are satisfied provided that v is small enough. 


Observe that (5.49) and (5.54) yield that the measure m is concentrated on [—1,1], The function 
b is therefore defined for m-a.e. zEK, even if the function T is only defined in ] — p , p[. 

Proof. We proceed according to the following steps. 


Q Step 1: we establish (3.30). We first point out that, by combining (5.10), (5.49) and (5.50), we 
can conclude that the following estimate holds for every iER: 


I Ufix) - Ufi = 


\fJ{x + gz) — Uj\fi(z)dz 


Ufix) — Uj <f>(z)dz 

Jr 

This implies that, if Uq E 53, then, since r < e 3 owing to (4.20c), we have 

\U 0 (x) - Ufi < | U 0 (x) - U,{x )| + | Ufix) — Ufi <r + 0( l)e 3 < e, 


< 0( l)e 3 


namely (3.30) holds true. 

0 Step 2: we establish the first inequality in (3.31a). We fix z E [—1,1] and we first point out 
that, owing to the explicit expression (4.18) of T, we have 

(5.55) b(z) > min{C c , C™} ■ « - <-i) = Cc • « - <-i) > 0 
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owing to (4.20). Also, for every z € [—1,1], we have that, owing to (5.13a), 

b(z) ^ (Xi Xi— \)£w — hvCw — & 

provided that v is sufficiently small because h v —> 0 + when v —> 0 + . This concludes the proof of 
the first inequality in (3.31a). 

Q Step 3: we establish the second inequality in (3.31a). We combine the explicit expression (4.17) 
of V with the inequality \V'\ < 0( 1)70 1 (see § 5.2) and we obtain that 


0 < £i(z) < 0(1)(x^ — Xj_ x )rj 1 < 0(1) i] 1 h u , for i = 1,2,3 and z E [—1,1] 

and this, jointly with ( 5.13a| ) and (5.55), implies the second inequality in (3.31a) provided that v 
is sufficiently small. 

Q Step 4: we establish ( 3.31d[ ). We first point out that, owing to ( |4.19| ) and ( 5.53[ ), 

U « + qz)-U «_! + sz) = [V « + qz) - V «_! + ?«)] + « + qz) - ® «_! + ?«)] 

= -b(z)fu - b(z)r 2I + b(z)r 3I 


+ d 3 


&(z),D 2 -U*),Di -£i(z),V~(z) 


-V~(z) 


By integrating the above equality with respect to the measure m and by recalling (3.31c) and (5.50) 
we conclude that 


(5.56) 


U q (x^) -U q (x IJ l _ l ) = —bf\j - br 2I + br 3I 


+ J r { D 3 ^(z),D 2 -£ 2 (z),D 1 -£i(z),v (z) -V (z)Jdm(^). 


Next, we recall (5.52) and we infer that 

(5.57) U+-U- = U 0 (xU)-Uo(x?) = U,(xU)-U^)-U,(xU) + Uo(xU)-Uo(x^) + U,(x^). 


Owing to ( |5.51[ ), if [To E 53, then 

Co«-i) - UM-d - 0 (xn - U«)] | < r(x\ - 
By recalling (5.57) with (5.51) and (|5.56[) we conclude that, if Uq E 35, then 


I U + — U + bru + bf 2I - br 3I — 


d 3 


£, 3 (z),D 2 -£ 2 (z),Di -£i( z),V (z) -V (z)\dm(z) 


< r« - - X0) < t 

To achieve the last two equalities we have used (4.20) and ( 5.55| ). This establishes (3.31d). 

0 Step 5: we establish (3.31b). We first recall that the first derivative of V satisfies \V'\ < 0(l)r] 
and we infer that the same bound holds for \U'\. By recalling (5.50), (5.49) and the inequality 
? < e 2 r/Cc we conclude that 

r 1 


-l 


\U(xU) -Us(xU)\ = 




[C/«_i) - U(x\_ 1 +qz)]<i>(z)dz 


< 


J 0 (l)rj 1 scf)(z)dz < 0(1 )e 2 ( c 
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By using again the inequality \V'\ < 0(l)r] , we infer that, for every Uq E IB and every z such 

that \z\ < 1 we have 


(5.58) 


I V~(z) -U~\<\V (xU + w) - E «_,) | + |E (xU) ~ U «-r) I 

+ |U (<_,) - % {xU) I + 1^ «-i) - (<-i) I 

<0(l)r/ _1 ?|2;| + 0(l)eC«,^ + 0(l)e 2 Cc + r < 0(l)e 2 r/. 

To establish the last inequality we have used (4.20), (5.14) and the inequality q < e 2 ry£ c . Next, we 


recall the explicit expression (4.17) of V and (5.53) and we conclude that 
(5.59) 0 < £1(2:) + £2 (z) +6(2:) < 0(1 )r]^ 1 C 1 (^ [x^ +sz,x" + sz\ n . 


By using the explicit expression (4.18) of T we infer that, for every 2 such that \z\ < 1, 


b(z) = Cc ( [<_ 1 + sz, xi + qz\ n 94 c ) + (w C} ( [x^ + <^, + <^] n 

> C™ ( [x-Li + x-' + n . 


By combining the above formula with (5.58) and (5.59) we eventually arrive at (3.31b). This 
concludes the proof of the lemma. □ 


5.8.2. Proof of Proposition 5.6: conclusion. In this paragraph we complete the proof of Proposi¬ 
tion |5.6[ which has Theorem |1.2| as a direct corollary. 

We consider the Baiti-Jenssen system (2.6) and the set 53 defined as in ( |5.51[ ). We have to show 
that conditions [ij) , |n|) and iii) in the statement of Theorem 1.2 are satisfied. Condition [I]) is satisfied 


owing to the considerations in § 2.2.1 Condition ii| is satisfied: indeed 


53 is nonempty since it contains (7 ? , because E <S(M). Also, 

53 is open in the topology of <S(M), which is stronger than the strong W 1,oa topology. 


We are left to show that condition |iii|) is also satisfied. The reason why condition condition iii|) is 
satisfied is because the proof Proposition 


4.4 


continues to work if we replace the function U with 
the function U q . provided that ? < e 2, r]Q. To see this, we first fix Uq E 53 and we introduce its 
wave front-tracking approximation by arguing as in § 5.3.1 Next, we point out that, owing to 


Lemma 5.7 we can apply Lemma 3.5 This implies that the same conclusions as at the end of 


5.3.2 and § 5.3.3 hold true. This in turn implies that all the analysis in § 5.4-§ 5.7 applies. We 


can infer that Proposition 4.4 holds true if we replace U with and hence we conclude the proof 


of Proposition 5.6 and Theorem 1.2 


Notation 

General mathematical symbols 

C N : the Lebesgue measure on 

0(1): any function satisfying 0 < c < 0(1) < C for suitable constants c, C > 0. The precise value 
of C and c can vary from line to line 

<S(M): the Schwartz space of rapidly decreasing functions, endowed with the standard topology (see 
for instance [23] p.133] for the precise definition) 

||-|liyioo : the standard norm in the Sobolev space W loa 

TotVart/: the total variation of the function U : M —> M. N , see [3, § 3.2] for the precise definition 
z\ ■ Z2 : the Euclidian scalar product between the vectors Z\, Z2 E 
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D- L [a, U\. the i- wave fan curve through U 

F': the hrst derivative of the differentiable function F : M —> M. N 
F(x±): the left and right limit of the function F at x (whenever they exist) 
Ri[s,U]: the z-rarefaction curve through U 
U]: the z-shock curve through U 
W 1 ’ 00 : the space of Lipschitz continuous functions 
a.e. (t,x): for £ 2 -almost every (t,x) 
a.e. x: for /^-almost every x 


See 


( 2 . 1 ) 


See 

See 


S|2T 


Symbols introduced in the present paper 

6 : a strictly positive parameter 

if. the perturbation parameter in the flux function F v 
Xi(U ): the z-th eigenvalue of the Jacobian matrix JF^ 
q, p: strictly positive parameters 
Jig,, 91 r : open subsets of M 

fi v \ the threshold for using the accurate Riemann solver 
v, h v : parameter and mesh size for the wave front-tracking approximation 
oj: a strictly positive parameter 
T: the function ^ : M —> M 3 


See fl4.16| ), ( |4.20aD 
See Q, jOPbl l 
See Q 
See (J4d6j), flOPa] ) 
See fl4J5| ) 
See 
See 8 


2.1 


5.3.1 


See (4.4), (4.20b) 


U', U". U* and U**: fixed states in M 3 

Vi', the strength of a shock i 

e: a strictly positive, sufficiently small parameter 

fi(U ): the z-th right eigenvector of the Jacobian matrix JF V 

T : the strictly positive interaction time 



See § 

$ 4.3.2, 4.3.3 

(4.20e) 

See § 4.3.3 

Remark 5.1 


See Page [ 6 ] 
See Proposition (|4.4[) 


U, Ur', the function U : 


and its mollification 


See 



See fl4.18D , ( |4.20c| ), ( [5.1aD 
See (|47l8|), (IOo3, d57lab 


See Lemma 4.1, Remark 5.1 


2 . 2.1 


( c : a strictly positive parameter 
( w : a strictly positive parameter 

q: a strictly positive parameter that we fix equal to 20 
r: a strictly positive parameter See (4.20c), (5.1), (4.21) 

u, w, v : the hrst, second and third component of the vector-valued function U See § 

U u : wave front-tracking approximation of the admissible solution U See 

Uq '. wave front-tracking approximation of the initial datum Uq 
V : the Lipschitz continuous function V : M —> M 3 
W: the piecewise constant function W : M —> M 3 
x ”: mesh points for the wave front-tracking approximation 


2.1 


See 8 2.1 


See ( |4.17D 
See p~T| ) 
See (5.13a) 


Ackowledgments. This work was originated by a question posed by Tai Ping Liu, whom the 
authors wish to thank. Part of this work was done when the hrst author was affiliated to the 
University of Oxford, which she thanks for the stimulating scientific environment. The second 
author thanks the University of Oxford for supporting her visits, during which part of this work was 
done. Both authors are members of the Gruppo Nazionale per l’Analisi Matematica, la Probability 
e le loro Applicazioni (GNAMPA) of the fstituto Nazionale di Alta Matematica (INdAM) and 



















































































































































46 


L. CARAVENNA AND L. V. SPINOLO 


are supported by the PRIN national project “Nonlinear Hyperbolic Partial Differential Equations, 

Dispersive and Transport Equations: theoretical and applicative aspects”. 

References 

[1] Adimurthi, Shyam Sundar Ghoshal, and G. D. Veerappa Gowda. Structure of entropy solutions to scalar con¬ 
servation laws with strictly convex flux. J. Hyperbolic Differ. Equ., 9(4):571—-611, 2012. 

[2] Luigi Ambrosio and Camillo De Lellis. A note on admissible solutions of ID scalar conservation laws and 2D 
Hamilton-Jacobi equations. J. Hyperbolic Differ. Equ., 1 (4):813—826, 2004. 

[3] Luigi Ambrosio, Nicola Fusco, and Diego Pallara. Functions of bounded variation and free discontinuity problems. 
Oxford Mathematical Monographs. The Clarendon Press, Oxford University Press, New York, 2000. 

[4] Fabio Ancona and Kliai T. Nguyen. In preparation. 

[5] Paolo Baiti and Helge Kristian Jenssen. Blowup in L°° for a class of genuinely nonlinear hyperbolic systems of 
conservation laws. Discrete Contin. Dynam. Systems, 7(4):837-853, 2001. 

[6] Stefano Bianchini and Alberto Bressan. Vanishing viscosity solutions of nonlinear hyperbolic systems. Ann. of 
Math. (2), 161(1):223—342, 2005. 

[7] Stefano Bianchini and Laura Caravenna. SBV regularity for genuinely nonlinear, strictly hyperbolic systems of 
conservation laws in one space dimension. Comm. Math. Phys., 313( 1): 1 -33, 2012. 

[8] Alberto Bressan. Hyperbolic systems of conservation laws. The one-dimensional Cauchy problem, volume 20 of 
Oxford Lecture Series in Mathematics and its Applications. Oxford University Press, Oxford, 2000. 

[9] Alberto Bressan and Rinaldo M. Colombo. Decay of positive waves in nonlinear systems of conservation laws. 
Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 26(1): 133-160, 1998. 

[10] Alberto Bressan and Philippe G. LeFloch. Structural stability and regularity of entropy solutions to hyperbolic 
systems of conservation laws. Indiana Univ. Math. J., 48(l):43-84, 1999. 

[11] Laura Caravenna and Laura V. Spinolo. New interaction estimates for the Baiti-Jenssen system. Preprint, 
available on arXiv. 

[12] Constantine M. Dafermos. Regularity and large time behaviour of solutions of a conservation law without 
convexity. Proc. Roy. Soc. Edinburgh Sect. A, 99(3-4):201-239, 1985. 

[13] Constantine M. Dafermos. Wave fans are special. Acta Math. Appl. Sin. Engl. Ser., 24(3):369-374, 2008. 

[14] Constantine M. Dafermos. Hyperbolic conservation laws in continuum physics, volume 325 of Grundlehren der 
Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin, 
third edition, 2010. 

[15] James Glirnm. Solutions in the large for nonlinear hyperbolic systems of equations. Comm. Pure Appl. Math., 
18:697-715, 1965. 

[16] James Glimm and Peter D. Lax. Decay of solutions of systems of nonlinear hyperbolic conservation laws. Memoirs 
of the American Mathematical Society, No. 101. American Mathematical Society, Providence, R.I., 1970. 

[17] Jonathan Goodman and Zhou Ping Xin. Viscous limits for piecewise smooth solutions to systems of conservation 
laws. Arch. Rational Mech. Anal, 121(3):235-265, 1992. 

[18] Helge Holden and Nils Henrik Risebro. Front tracking for hyperbolic conservation laws, volume 152 of Applied 
Mathematical Sciences. Springer-Verlag, New York, 2002. 

[19] S. N. Kruzkov. First order quasilinear equations with several independent variables. Mat. Sb. (N.S.), 81 
(123):228-255, 1970. 

[20] P. D. Lax. Hyperbolic systems of conservation laws. II. Comm. Pure Appl. Math., 10:537-566, 1957. 

[21] Tai Ping Liu. Decay to N- waves of solutions of general systems of nonlinear hyperbolic conservation laws. 
Comm. Pure Appl. Math., 30(5):586-611, 1977. 

[22] O. A. Oleinik. Discontinuous solutions of non-linear differential equations. Uspehi Mat. Nauk (N.S.), 12(3(75)):3- 
73, 1957. 

[23] Michael Reed and Barry Simon. Methods of modern mathematical physics. I. Functional analysis. Academic 
Press, New York-London, 1972. 

[24] David G. Schaeffer. A regularity theorem for conservation laws. Advances in Math., 11:368-386, 1973. 

[25] Denis Serre. Systems of conservation laws. 1. Cambridge University Press, Cambridge, 1999. Hyperbolicity, 
entropies, shock waves, Translated from the 1996 French original by I. N. Sneddon. 



A COUNTER-EXAMPLE TO THE REGULARITY OF SYSTEMS 


47 


[26] Eitan Tadmor and Tamil' Tassa. On the piecewise smoothness of entropy solutions to scalar conservation laws. 
Comm. Partial Differential Equations, 18(9-10):1631 1652, 1993. 

[27] Tao Tang, Jinghua Wang, and Yinchuan Zhao. On the piecewise smoothness of entropy solutions to scalar 
conservation laws for a larger class of initial data. J. Hyperbolic Differ. Equ., 4(3):369-389, 2007. 

L.C. Dipartimento di Matematica, Universita degli Studi di Padova, Via Trieste 63, 35121 Padova, 
Italy 

E-mail address', caravenna@math.unipd.it 

L.V.S. IMATI-CNR, via Ferrata 1, 1-27100 Pavia, Italy. 

E-mail address', spinolo@imati.cnr.it 



